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Abstract 


The Relativistic Continuum Hartree-Bogoliubov (RCHB) theory, which properly takes into account the 
pairing correlation and the coupling to (discretized) continuum via Bogoliubov transformation in a micro- 
scopic and self-consistent way, has been reviewed together with its new interpretation of the halo phenomena 
observed in light nuclei as the scattering of particle pairs into the continuum, the prediction of the exotic 
phenomena — giant halos in nuclei near neutron drip line, the reproduction of interaction cross sections and 
charge-changing cross sections in light exotic nuclei in combination with the Glauber theory, better restora- 
tion of pseudospin symmetry in exotic nuclei, predictions of exotic phenomena in hyper nuclei, and new 
magic numbers in superheavy nuclei, etc. Recent investigations on new effective interactions, the density 
dependence of the interaction strengthes, the RMF theory on the Woods-Saxon basis, the single particle 
resonant states, and the resonant BCS (rBCS) method for the pairing correlation, etc. are also presented in 
some details. 
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1 Introduction 


Currently nuclear physics is undergoing a renaissance as evidenced by the fact that worldwide there are lots 
of Radioactive Ion Beam (RIB) [1] facilities operating, being upgraded, under construction or planning to be 
constructed, e.g., the Cooling Storage Ring (CSR) project at HIRFL in China will be completed in 2005 [2], 
the RIB Factory at RIKEN in Japan will begin to operate at the end of 2006 [3], the FAIR project at GSI in 
Germany was approved in 2003 [4], and the construction of the Rare Isotope Accelerator (RIA) is considered 
to be of the highest priority of all physics disciplines in the US [5], etc. These new facilities together with 
developments in the detection techniques have changed the nuclear physics scenario. It has come into reality 
to produce and study the nuclei far away from the stability line — so called “EXOTIC NUCLEI” [6-12]. 


New exciting discoveries have been made by exploring hitherto inaccessible regions in the nuclear chart |6— 
12]. For examples, a new phenomenon, halo — a state in which nucleons spread like a thin mist around the 
nucleus, was first discovered in HLi with RIB in 1985 [13] and later in many other light exotic nuclei as well; 
exotic nuclei exhibit other interesting phenomena such as the disappearance of traditional shell gaps and the 
occurrence of new ones, which result in new magic numbers [14]; etc. The exotic nuclei also play important 
roles in nuclear astrophysics as well, as their properties are crucial to understand stellar nucleosynthesis. 

The current nuclear models are mainly based on the knowledge obtained from the nuclei near the 6- 
stability line. For instances, the cornerstones for the edifice of modern nuclear physics include shell model [15] 
and collective model [16,17], which are respectively based on the magic numbers (the stability of some nuclei 
compared with their neighbors) and the incompressibility of the nuclear matter. Therefore the change of 
magic numbers and the unprecedented low density nuclear matter in halo nuclei have shaken the foundation 
of nuclear physics. New innovative nuclear models are needed to describe the exotic nuclei characterized by 
weakly binding and low density. 

The relativistic mean field (RMF) theory [18] has received wide attention due to its successful description 
of lots of nuclear phenomena during the past years [19-23]. In the framework of the RMF theory, the nucleons 
interact via the exchanges of mesons and photons. The representations with large scalar and vector fields in 
nuclei, of order of a few hundred MeV, provide simpler and more efficient descriptions than nonrelativistic 
approaches that hide these scales. The dominant evidence is the spin-orbit splittings. Other evidence includes 
the density dependence of optical potential, the observation of approximate pseudospin symmetry, correlated 
two-pion exchange strength, QCD sum rules, and more [24]. The relativistic Bruekner-Hartree-Fock theory 


H ` and RMF theory with density-dependent coupling constants extracted from it can reproduce well the nuclear 
ON saturation properties (the Coester line) in nuclear matter [25,26]. Furthermore the RMF theory can reproduce 
— better the measurements of the isotopic shifts in the Pb-region [27], give more naturally the spin-orbit potential, 
c ` the origin of the pseudospin symmetry [28,29] as a relativistic symmetry [30-32] and spin symmetry in the 
JO anti-nucleon spectrum [33], and be more reliable for nuclei far away from the line of 8-stability, etc. Obviously, 
— RMF is one of the best candidates for the description of exotic nuclei. 


In order to describe exotic nuclei, the pairing correlation and the coupling to continuum, which are extremely 


=) crucial for the description of drip line nuclei, must be taken into account properly [34]. The continuum effect 


N is commonly taken into account in the Hartree-Fock-Bogoliubov (HFB) [35] or relativistic-Hartree-Bogoliubov 


e (RHB) [21] approaches. In most of these calculations the continuum is replaced by a set of positive energy states 
= determined by solving the HFB or RHB equations in coordinate space and with box boundary conditions |36- 
“ 38]. Recently the HFB equations were also solved with exact boundary conditions for the continuum spectrum, 
- both for a zero range [39] and a finite range pairing forces [40]. These resonant continuum HFB or RHB 
== results with exact boundary conditions are generally close to those with box boundary conditions [39,41]. The 
“= extension of the RMF theory to take into account both bound states and (discretized) continuum via Bogoliubov 


" transformation in a microscopic and self-consistent way, i.e., the Relativistic Continuum Hartree-Bogoliubov 
(RCHB) theory has been done in Refs. [36,42]. The RCHB theory was very successful in describing the ground 
state properties of nuclei both near and far from the f-stability line. By using a density-dependent zero 
range interaction, the halo in "Li has been successfully reproduced in this self-consistent picture. Remarkable 
successes of the RCHB theory include the new interpretation of the halo in HLi [36], the prediction of the 
exotic phenomena as giant halos in Zr (A > 122) [43] and Ca (A > 60) [44] isotopes, the reproduction of 
interaction cross section and charge-changing cross sections in light exotic nuclei in combination with the 
Glauber theory [45,46], better restoration of pseudospin symmetry in exotic nuclei [31,32], and predictions of 
exotic phenomena in hyper nuclei [47] and new magic numbers in superheavy nuclei [48], etc. 

'The main purpose of the present manuscript is to review the RCHB theory and its applications for exotic 
and superheavy nuclei in spherical case. Some other related topics will also be covered briefly. In Section 2, 
the formalism, the numerical solutions and the effective interactions for the RMF theory are presented. In 
Section 3, the pairing correlations and the approaches for pairing are sketched. Discussion on the continuum 
states follows where the resonant BCS method and several methods to obtain single particle resonant states 
are briefly reviewed. Then the detailed formalism of the RCHB theory together with the discussion on the 
effective pairing interactions are given. In Section 4, the applications of the RCHB theory to the properties 
of exotic nuclei are presented, such as the binding energies, particle separation energies, the radii and cross 


sections, the single particle levels, shell structure, the restoration of the pseudo-spin symmetry, the halo and 
giant halo, and halos in hyper nuclei, etc. In Section 5, the predictions of new magic numbers in superheavy 
nuclei are presented. Finally a brief summary and perspectives are given in Section 6. 


2 Relativistic mean field theory 


In this section, we present the formalism, the numerical solutions and the effective interactions for the RMF 
theory as well as its application for nuclear matter. In the first subsection, the effective Lagrangian density and 
equations of motion for the nucleon and the mesons are given. The numerical solutions of the Dirac equation 
for finite nuclei are discussed in the second subsection. Subsequently, the effective interactions with nonlinear 
self-coupling meson fields and density dependent meson-nucleon couplings, and their influences on properties 
of nuclear matter are discussed. 


2.1 The general formalism 


The RMF theory describes the nucleus as a system of Dirac nucleons which interact in a relativistic covariant 
manner via meson fields. The meson fields are treated as classical fields. In the simplest RMF version, i.e., 
the o-w model [18], the mesons do not interact among themselves, which leads to too large incompressibility 
in nuclear matter. Therefore a nonlinear self-coupling of the o-field was proposed [49]. In order to repro- 
duce the density dependence of the vector and scalar potentials of the Dirac-Brueckner calculations [26], the 


~ nonlinear self-coupling of the w-meson was found to be necessary [50]. Recently the nonlinear self-coupling of 
Kei the isovector p-meson was also introduced to improve the density-dependence of the isospin-dependent part of 
e the potentials [51]. Within the present scheme, the isoscalar-scalar o-meson provides the mid- and long-range 


) attractive part of the nuclear interaction whereas the short-range repulsive part is provided by the isoscalar- 
vector w-meson. The photon field A“(a) accounts for the Coulomb interaction while the isospin dependence 


=) of the nuclear force is described by the isovector-vector p-meson. The z-meson field is not included because it 
| does not contribute at the Hartree level. In principle, other mesons apart from c, w, and p may also contribute 
™ to the nuclear interaction as well, e.g., the isovector scalar ö-meson, which was suggested in Ref. [52] and also 
= ` on the basis of a relativistic Brueckner theory in Refs [53-56], has been used in some structure calculations for 
^ = exotic nuclei [55], nuclear matter [57] and stellar matter [58]. However, as the RMF theory is only an effective 
‚= theory therefore it is expected that the contributions from other mesons can be effectively taken into account 
Sc by adjusting the model parameters to the properties of nuclear matter and finite nuclei. With various versions 


(© of the nonlinear self-couplings of meson fields, the RMF theory has been used to describe lots of nuclear phe- 
` = nomena during the past years with great successes [19-21]. In order to avoid the problem of instability for 


€^ the nonlinear interactions at high densities, the RMF theories with density dependence in the couplings are 
) developed as well [26, 51, 59-62]. 
The Lagrangian density of the RMF theory can be written as 


—-|. =f Ze 1-73 
L= v iy"0,— M — goo — gui" wy — go Y" T - Pu — ey" Ay Y 
lou lou lw. g ~ (1) 
+ 59 oô a — U,(o) — 72 Ou + Ugly) — 7E iau Up Pu) 
1 
u PP" Fw, 


where M and m;(g;) (i = 0,w, p) in the following are the masses (coupling constants) of the nucleon and the 
mesons respectively and 


Q^" = Of — wh, (2a) 
RY = Oe Or, (2b) 
p" = ah AY _ "AM (2c) 


are the field tensors of the vector mesons and the electromagnetic field. We adopt the arrows to indicate vectors 
in isospin space and bold types for the space vectors. Greek indices u and v run over 0, 1, 2, 3 or t, z, y, z, 
while Roman indices i, j, etc. denote the spatial components. 


The nonlinear self-coupling terms U5(c), U,(w,), and Up(p,,) for the o-meson, w-meson, and p-meson in 
the Lagrangian density (1) respectively have the following forms: 


. Pow mu 1 2.1. 4 
Ls = SEN + 3920 + 7930 ; (3a) 
1 1 
U,lay,) = g niu y + 70s (Wwy), (3b) 
4 EE PET 
Dë = mpg": Put dal" Au”. (3c) 


From the Lagrangian density (1), the Hamiltonian operator can be obtained by the general Legendre 


H = fez H= [ez 5 ni (a) Oe) — £(x)| , (4) 
where the conjugate momenta of the field operators ¢; (Qi = V, 0, wv, py, Av) are defined as 
(x) = PUN (5) 
'Then the Hamiltonian density of the system can be easily obtained as 
H —H[v] + Hlo] + ev] + Ho] + [A] 
= [-i7'0; + M] b+ goo) + gown Y" b + gp : Vy" T + eA” 


1 bL 1 1.; 

+ 50700 = 39 cic + U, (c) = 7000, + nit — TC 
ts - lo, = 1 1; 

= read Rov + LET - Riv — Dell" qe Fo + PF. 


The relativistic mean field theory is formulated on the basis of the effective Lagrangian (1) with the mean 


transformation 


1-73 


Y 


Y - field approximation, i.e. the meson fields are treated as classical c-numbers. With this approximation, all the 
— quantum fluctuations of the meson fields are removed and the nucleons are described as independent particles 
^. moving in the effective meson and photon fields. Therefore the nucleon field operator can be expanded on a 


> complete set of single-particle states as, 


da) =) dala) (7) 


‘= where c, is the annihilation operator for a nucleon in the state a of the Dirac fields and v, is the corresponding 
_@ single-particle spinor. The operator c, and its conjugate ch satisfy the anticommutation rules 


(ea, cl) = ab and edu Il. di e. (8) 
Confined to the single-particle states i with positive energies, i.e. the no-sea approximation, the ground state 
of the nucleus can be constructed as, 


A 
|j -[[d]o with (8/6) =1, (9) 
i=1 


where |0) is physical vacuum. 
With the ground state (9) and the mean field approximation, the energy functional, i.e. the expectation 
value for the Hamiltonian (6) is obtained as 


Ermr (P, $) =(®|H|®) 


1 
= / d’xTr l (7 PEM + Goo + gue yu + gef Bu" + sell: na) Aya) d 


1 1 1 = E 1 
+ / dix [- gods + Ug (e) + =O!” — Us wy) + FRM -Ru — U (Ru) + iP! Fa} 
(10) 


where ¢ = 0, wy, Py, A, and the density matrix p is defined as 


pg = (@lWlelvje;|®) = vivió; (11) 


For system with time reversal symmetry, the space-like components of the vector fields vanish. Furthermore 
one can assume that in all nuclear applications the nucleon single-particle states do not mix isospin, i.e. the 
single-particle states are eigenstates of 73, therefore only the third component of p, survives. Stationarity 
implies that the nucleon single-particle wave function can be written as 


dile) = e" (x). (12) 
where e; is the single-particle energy. Accordingly the density matrix is reduced as, 
pij = v (x! )bi(x) big. (13) 


Altogether there remain only the meson fields c, wo, PÈ, Ao which are time independent. 

The equations of motion for the nucleon and the mesons can be obtained by requiring that the energy 
functional (10) be stationary with respect to the variations of p and ¢. More explicitly, the stationary condition 
reads 

ô (Ermr(p, $) — Tr(£p)) = 0, (14) 
_ where € is a diagonal matrix, whose diagonal elements are the single particle energies e;(i = 1,--- , N) intro- 
La duced in Eq.(12), and N is the number of eigenstates. 
N _ Using the variation dp with respect to v;, the stationary condition (14) leads to the Dirac equation 


(a-p+B(M +S) +V) vi(x) = eivi(x) (15) 


== for the nucleon and the Klein-Gordon equations for sigma, omega, rho, and the photon: 


eed m) = deis (16a) 
—V?ug + U, (wo) = Jup, (16b) 
-V° +U) = gem (16c) 
MI An = epe. (16d) 
= The scalar potential S and vector potential V in equation (15) are respectively: 
S(x) = 9o0(x), (17a) 
E V(x) = guteo(x) + gpT3P (x) + 2e (1 — n3) Aol). (17b) 


= While the scalar density (ps), the baryon density (p,), the isovector density (p3), and the charge density (pe) 
© in the Klein-Gordon equations (16a-16d) are respectively, 


ps = Trip], (18a) 
p. = Tele], (18b) 
ps = Tr[rsp], (18c) 
1 
pe = VT - nl. (18d) 
The total energy of the system can be obtained from the energy functional (10) as, 
1 1 
E= Io, {m [(a -p + 8M) p] + 5T: [9.800 + gutoop + gorapop + zei — SE } 
(19) 


» f dx [Us (o) — Us (wo) — Heft - 5 / da [oU (0) — uU! (wo) — gtt. 


In the density dependent RMF approach, where the nonlinear self-couplings for the c, w, and p mesons in the 
Lagrangian density are respectively replaced by the density dependence of the coupling constants gc (p), gu (p), 
and gj (p), an additional term, i.e. the rearrangement term, will appear in the Dirac equation (15) [26,51,59-62]. 


2.2 Numerical algorithm for spherical nuclei 


The harmonic oscillator basis has served as a very useful tool in nuclear structure study. Normally, the equations 
of motion for nucleons moving in a mean field are solved by expanding them on the harmonic oscillator (HO) 
basis [15-17, 35,63]. However, for exotic nuclei with large spacial extension, e.g., halo nuclei, it is not justified 
to work in the conventional harmonic oscillator basis due to its localization [36,64—66]. Instead, one can choose 
to work either in the coordinate space, or improve the asymptotic behavior of the HO wave function, or adopt 
other basises which have a correct asymptotic behavior, for example, the Woods-Saxon basis. 

In this subsection, we will focus on the numerical solution of the RMF for spherical nuclei. Due to the 
special spacial symmetry, both the Dirac equation for the nucleon and the Klein-Gordon equations for the 
mesons and the photon become radially dependent only, thus facilitating much the solution of the coupled 
equations. The formalism for the spherical relativistic Hartree (SRH) theory will be briefly presented. We 
then review the SRH theory in different basis, including Finite Element Method (FEM) in coordinate space, 
transformed harmonic oscillator basis, and the Woods-Saxon basis. The application of the SRH theory to 
doubly magic nuclei follows. 


2.2.1 Spherical relativistic Hartree theory 


Starting from the Eqs. (15) and (16a-16d) given in the previous subsection, one derives the coupled radial 
equations for spherical nuclei, i.e., the radial Dirac equation and radial Klein-Gordon equations. 
For spherical nuclei, the Dirac spinor which is the expansion coefficient Ya(r) (a = (o, &, m)) in Eq. (12) 


— (the coordinate from x has been changed to r to reflect the spherical symmetry) is characterized by the angular 


«i momentum quantum numbers «(l,j), m, the parity, the isospin t = +1/2 (“+” for neutrons and “—” for 


= d protons) and the radial quantum number o and has the form 


Gale) 


r 


Vom E t) = F*(r) 


GC P)Xta (t), (20) 


C * with G5 (r)/r and Fz(r)/r the radial wave functions for the upper and lower components and Val, $) the 
spinor spherical harmonics [67]. Substituting Eq. (20) into the Dirac equation (15), one can deduce the radial 


_ > Dirac equations as 


Cat? = (-z + d FS TM + S(r)+V(r)) G$, (21a) 
GE = (+5 + 3 Gs - (M + S(r) -V(r) Fe, (21b) 
with the scalar and vector potentials 
S(r) = Hot, (22a) 
1 
Vir) = guwo + IpT3P + zei — 73) Ao. (22b) 


The meson field equations (16a-16d) simply become radial Laplace equations of the form 


d^ 2d 
(3-25 + md) an sun. (x) 
where mg are the meson masses for ® = c, c, p and zero for the photon. The source terms are 


—geps(r) — gaa?(r) — gac?(r), for the o—field, 
Salt) = JuPu(r) — eawi(r); for the w—field, 

gopa(r) — daleg(r)P?, for the p—field, 

epc(r), for the Coulomb field, 


(24) 


with 


A 
er) = en Gf - OP), (25a) 
i P 
e) = qs RR, (25) 
: Ni 
pl) = geo HRO, (250) 
1 &1 , ; 
per) = wén IFO IG? + ROP). (254) 


The procedures of solving these coupled equations are as the following: a) with a set of estimated meson 
and photon fields, the scalar and vector potentials in Eqs.(22) are calculated and the radial Dirac equation 
is solved; b) the so obtained nucleon wave functions are used to give the source term of each radial Laplace 
equation for mesons and the photon; and c) the new meson and photon fields obtained from solving these 
Laplace equations will be used to replace the fields in step a). This procedure is iterated until a demanded 
accuracy is achieved. 

With the above procedures, there are two methods to solve the coupled equations (21-25). One is in 


=> coordinate space with the shooting method [68], i.e. spherical relativistic Hartree theory in r space (SRHR), 
, and the finite element method [42, 69-71] (SRHFEM). Another is in configuration space, e.g., the harmonic 
* oscillator basis [72] (SRHHO), the transformed harmonic oscillator basis [64] (SRHTHO) and the Woods-Saxon 
„m basis [73] (SRHWS). The readers are referred to Ref. [68] and Ref. [72] for SRHR and SRHHO, respectively. 
: > In the following the numerical techniques for the solution of the Dirac equation in the finite element method, 
= the transformed harmonic oscillator basis and Woods-Saxon basis will be briefly introduced. 


c2 2.2.2 Finite element method 


— A convenient procedure for the coordinate space discretization of the Dirac equations (21) is provided by Finite 
—— Element Method (FEM) [42,69-71]. Similar as in the r space [68], the Dirac equations (21) are solved in a box 


æ With a box size R and proper boundary condition. The radial wave functions G% and FẸ are discretized at 
= N +1 points, i.e., rj = 0, r2 = Ar, ... , rN41 = N - Ar, where Ar = R/N. The wave functions thus obtained 


= become tabulated values: GA (rg),G (r1), ... , G&(rw41) and Fira), Ein), ... , FE(r N41). 


The idea of FEM is as following: if the Ar — R/N is small enough, or equivalently N is large enough, the 


' wave function G^ (r) in the element r € [r;,r;+1] can be well approximated by G(r;) and G$ (ri+1) together 


with some simple analytic function, e.g., the linear, quadratic, cubic, or 4th order shape function [69]. 
Taking the linear shape function, 


a = 1- ` 
nalr) p (26) 
mr) — p 
with p — = "i andre [ri, ri 41] as an example, the wave function G$ (r) can be well approximated by, 
T 
Gar) = Ga (ri)m(r) + GS (rii)m(r), rE [ri ri]. (27) 


Similarly the expression for F(r) can also be obtained. For r € [r;,r;+1], the Dirac equations (21) can now be 
written as: 


M S(r)*V(r)-e - +7 | GR(ri)na(r) + G&(risa)m(r) E (28) 
Zee M—S(r) + V(r) ea | \ Falranalr) + Eër) 


Ti+1 


Multiplying (G2 (ri)na(r)+G$ (riga)no(r) FS (ri)na(r) + FS (rizi)m(r)) from left and integrating with / r?dr, 


ri 
algebra equations for G$ (ri), Girl and FZ(r;i), F&(r;+1) are obtained. Repeating the same procedure for 
all the N elements, the Dirac equations (21) are transformed into a generalized eigenvalue problem. This 
generalized eigenvalue problem can be solved by the standard diagonaliztion algorithms and all the energies 
and wave functions discretized at r = rı,ra,...,rnı+ı can be obtained. It has been shown that FEM can 
provide very accurate solutions for the relativistic eigenvalue problem in the self-consistent mean-field approx- 
imation [42, 69-71]. 

Comparing with the shooting method in the r space [68], the FEM has the advantage that one can get 
all the energy e, and wave functions G% and F by a single diagonalization and it is straight forward to be 
generalized to the cases with nonlocal interaction in the pairing channel. However, in order to obtain the same 
accuracy as the shooting method, a huge matrix has to be constructed and it becomes more time consuming. 
One can also replace the linear shape function er) and mir) in Eq.(26) by the quadratic, cubic, or 4th order 
shape function, the procedures are the same and in principle the same accuracy can be achieved by the linear 
shape function with larger N. 


2.2.3 Transformed harmonic oscillator basis 


In order to modify the asymptotic behavior of the harmonic oscillator wave function at large r, the local- 
a scaling transformation method [74-76] has been introduced to construct the so called local-scaling transformed 
Z harmonic oscillator basis (THO) in Refs. [64,65]. 

A local-scaling point coordinate transformation (LST) is defined as 


r' =f(r) 2 £f(r). (29) 


yo, The transformed radius vector has the same direction f = r/|r|, while its magnitude r’ = f(r) depends on the 
=) scalar LST function. f(r) is assumed to be an increasing function of r, and f(0)=0. The corresponding LST 


wave function can be expressed as 


vr), Fra), Fra), (30) 


V (r1, r2, TA) = | 


i21 í 


> where V(r;,ro,...,rA) is an A-particle wave function normalized to unity. The local one-body density corre- 
“© sponding to an A-body wave function V is 


p(r) =A [ Wera. raros tra. (31) 


There exists a simple relation between the local density pr(r) associated with the LST function V ;, and the 
density p(r) which corresponds to the prototypical model function dr: 


2 r r 
pra) = FPL A en, (32) 


r2 


When the form of the density p(r) is known, Eq. (32) becomes a first order nonlinear differential equation for 
the LST function f and can be solved easily. For a system with spherical symmetry, pf, p, and f depend only 
on r=|r|, and Eq. (32) can be reduced to a nonlinear algebraic equation 


r f(r) 
at = f p(u)u?du. (33) 


The solution can be found subject to the boundary condition f(0) — 0. 
For shell-model or mean-field applications, one has to consider the case when the model many-body wave 
function is a Slater determinant 


: NE: 


V (ri, r2, ...,T A) 7a; 9e ler). (34) 
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Figure 1: Self-consistent RHB proton and neutron densities for the ground-states of "Ne and Ne. Density 
profiles calculated in THO are compared with FEM. Taken from Stoitsov et al. [64]. 


The single-particle functions ¢;(r) form a complete set. The LST wave function is defined by the transformation 
(30), and is written as a product state 


V r(vi, T2, ..., rA) VAI det lé; (rj). (35) 
ə of the transformed basis state T 
2 ð 
niet - LES ` Amen (36) 


In Refs. [64,65], the transformed harmonic oscillator basis (THO) has been derived by a local scaling-point 
transformation of the spherical harmonic-oscillator radial wave functions. The unitary scaling transformation 
produces a basis with improved asymptotic properties. The THO basis is employed in the solution of the 
relativistic Hartree-Bogoliubov (RHB) equations in configurational space [64]. As shown in Fig. 1, an expansion 
- of nucleon spinors and mean-field potentials in the THO basis reproduces the asymptotic properties of neutron 
densities calculated by FEM in the coordinate space. 


2.2.4 Woods-Saxon basis 
Woods-Saxon basis from the Schródinger equation (SWS basis) For the Schródinger equation with 


a spherical Woods-Saxon potential 


EM MET 
Vws(r) = 1+ e(r — Ro)/ao ’ maxs T 


oo, r > Rmax, 


where Rmax is introduced for practical reasons to define the box boundary. The eigenfunction can be written 
as nim (r) = Rni(r)Yim, (0,9) and its radial Schrödinger equation is derived as 


1 1d 4,d  LI(L+1) 
ru — —T — P À 
2M \r?dr dr r? 


) + Vest Rar) = ni yir). (38) 


Equation (38) is solved on a discretized radial mesh with a mesh size Ar. Rmax (Ar) is chosen large (small) 
enough to make sure that the final results do not depend on it. The radial wave functions thus obtained form 
a complete basis, 


Urine 1-1 101 s (39) 
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in terms of which the radial part of the upper and the lower components of the Dirac spinor in Eq. (21) are 
expanded respectively as 


Nmax 


s(r) 2 -i Kë Ut Fir) 


(40) 


The radial Dirac equation (21) is transformed into the WS basis as 


Lë (m ue) (A1) 
Cun Dm Jan Jan 


where the matrix elements are calculated as follows 


Rmax 
Amn = rdr Rmı(r) (V(r) + S(r) + M) Rni(r) 
Rmax = 
Bin = / r dr Rmı(r) (+4 u = J (r), 
0 dr r n (42) 
ER m 
m= T aripi" E " nl (T), 
Rimax 
Das = / rdrR.;(r) (V(r) - S(r) - M) Rr) 
0 


In practical calculations, an energy cutoff Ec (relative to the nucleon mass M) is used to determine the 


` cutoff of the radial quantum number nmax [73]. 


L Woods-Saxon basis from the Dirac equation (DWS basis) The radial Dirac equation (21) may be 
= solved in r space [68] with Woods-Saxon-like potentials for Vo(r) + So(r) [77] within a spherical box of the size 
ox, Rmax, together with the spherical spinor which gives a complete WS basis 


es Lanes 5, t) Ire ? Geen > So}, (43) 


e ` with n = 0,1,---, & = £1,42,---, and m = —jx,-++ e, ant, 5,t) takes the form of Eq. (20). In such 
` cases, states both in the Fermi sea and in the Dirac sea should be included in the basis for completeness. The 
` nucleon wave function (20) can be expanded in terms of this set of basis as 


Nmax 


Vom, =}. Can ess (5, 8, t), (44) 
n=0 
where nmax = nj + Nmax + 1 and the summation runs over Boe energy levels in the Fermi sea for 


0 € n € ntas and over negative energy levels in the Dirac sea for nf y +1 X n < Nmax. The negative energy 
states are obtained with the same method as the positive energy ones. In this WS basis, the Dirac equation 
(21) turns out to be 


Cane. SR > Bass = €alCam, M = 1, ttt e Nmax; (45) 
n=0 
with 
Hy, = Qn) [AV (E) + BAS(r)] lite) 


= n = drG? (r) [AV (r) + AS(r)] G(r) 
0 

4 n C" IFO (r) [AV (r) = AS(r)] E»), (46) 
0 
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Table 1: Effects of negative energy levels on bulk properties in SRHDWS for 160. The effective interaction 
for the Lagrangian is NLSH, Rmax = 20 fm, Ar = 0.1 fm, and Nj, = 25. For the initial Woods-Saxon like 
potentials, parameters in Ref. [77] are used except for Vo which is specified in the table. The left value in each 
entry gives the result without negative energy levels included and the right one that with NA = 25. Energy 


is in MeV and radius in fm. Taken form Ref. [73]. 


ax 


8.547 | 8.013 | 8.117 | 8.015 | 8.427 | 8. 
2.385 | 2.568 | 2.531 | 2.567 | 2.610 | 2.567 


where AV (r) = V(r) — Vo(r), AS(r) = S(r) — So(r) and the angular, spin, and isospin quantum numbers are 
omitted for brevity. 

In the expansion (44) of the nucleon wave function in the SRHDWS theory, one has to take into account 
not only the states in the Fermi sea but also those in the Dirac sea because these states form a complete basis 
together. The contribution from negative energy states for 16O is given in Table 1. It is found that, without 
including the negative energy levels, the calculated results will depend on the potentials for the basis. 

The contribution of negative energy states in the Dirac sea to the wave function can be calculated by 


N S, lc|? in the expansion (44) and it is around 10~4~~° (Note that the nucleon wave function is normalized 


) to one). However, such a small component from negative energy states in the wave functions contributes to 


( ` the physical observables such as E/A and rims by magnitudes of 1~10 % as can be seen from Table 1. 


. Comparisons between r space, harmonic oscillator basis and Woods-Saxon basis It is found that 
for stable nuclei, the SRHR, SRHSWS, SRHDWS, and SRHHO approaches are all valid and results from them 
) are in excellent agreement with each other [73]. However, for unstable nuclei near the neutron drip line, these 


( - methods differ from each other to some extent. 


In Fig. 2, the neutron density distribution of Ca from different SRH approaches are compared. With 

«^ the same box size, the density distribution from SRHR are almost identical with those from SRHWS, which 

‘@@ indicates the equivalence between SRHWS and SRHR. For brevity, only pn(r) from SRHR with Rmax = 35 fm 

is displayed which covers the curve corresponding to pn(r) from SRHWS with Rmax = 35 fm in Fig. 2. On the 

— other hand, p,(r) from SRHHO even with Nmax = 43 fails to reproduce the result of SRHR due to the well 
; known localization property of HO potential [64]. 

These results indicate that even the long tail (or halo) behavior in neutron density distribution for nuclei 
near the drip line can be reproduced quite well by the expansion method in the Woods-Saxon basis in the 
Schródinger framework (SRHSWS) and that in the Dirac framework (SRHDWS). This can be seen that the 
neutron density distribution obtained in SRHR is reproduced well in both SRHSWS and SRHDWS when 
enough box size is taken. 


2.3 Effective interactions 


In the Lagrangian density (1), there are meson masses Mo, Mw, Mp and meson-nucleon coupling constants 
Jo, Jw, gp together with the nonlinear self-couplings of the meson fields left to be determined. They are the 
nucleon-nucleon interactions in the RMF theory and in principle should be determined either by more funda- 
mental theories or by experiments. However, as the relativistic mean field theory is formulated on the basis of 
the above effective Lagrangian in connection with the mean-field and the no-sea approximations, it is difficult 
to determine these interactions microscopically. Instead, the masses and coupling strengthes of the mesons and 
the nonlinear self-couplings of the meson fields are determined by reproducing the properties of nuclear matter 
and a few doubly magic nuclei. Namely, they are effective interactions in the similar sense as their conventional 
counterparts. The effective interactions in the RMF theory could be determined by minimizing the least square 
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æ Figure 2: Comparison of density distributions for 72Ca from SRHR, SRHWS and SRHHO with NLSH. Ar = 
` 0.1 fm for SRHR and SRHWS. E = 75 MeV and Rmax = 20 (thick dot-dashed curve), 25 (thick long-dashed 


= ; curve), 30 (thick dashed curve) and 35 fm (thick dotted curve) for SRHWS. These sets of cutoff's correspond 
= to cutoff's in principal quantum number Nmax = 25 (thin dot-dashed curve), 31 (thin long-dashed curve), 37 


e (thin dashed curve) and 43 (thin dotted curve) which are used in SRHHO calculations. Taken form Ref. [73]. 


error as follows, 


N | exp T 


2525 i vena 


= 
i=l KR 


where a is the ensemble of the meson masses Mo, Mw, m, and the meson-nucleon coupling constants go, gw, gp 


“> together with the nonlinear self-couplings of the meson fields to be fitted, and y; " and c; are the experimental 
æ observables and corresponding weights. In general, the masses and charge radii of spherical nuclei near the 


^ B-stability line are adopted as observables in the least-square fitting procedure. 

-— Among the existing effective interactions for the RMF theory, the most frequently used are NL1 [78], 
>= NLSH [79], TM1 [50] and NL3 [80] with nonlinear self-couplings of mesons. Along the ß-stability line NL1 
gives excellent results for binding energies and charge radii, in addition it provides an excellent description 
of the superdeformed bands [81,82]. However, when moving away from the stability line the results are less 
satisfactory. This can be partly attributed to the large asymmetry energy J ~ 44 MeV. Moreover, the neutron 
skin thicknesses calculated with NL1 show systematic deviations from the experimental data [83]. 

In NLSH this problem was treated in a better way and the improved isovector properties have been obtained 
with an asymmetry energy of J ~ 36 MeV. Furthermore, NLSH seems to describe the deformation properties 
in a better way than NL1 does. However, NLSH produces a slight over-binding along the line of 6-stability 
and in addition it fails to reproduce successfully the superdeformed minima in Hg-isotopes in constrained 
calculations for the energy landscape. A remarkable difference between the two effective interactions are the 
quite different values predicted for the nuclear matter incompressibility [84], i.e., K = 212 MeV for NL1 while 
K = 355 MeV for NLSH [85,86]. As an improvement, NL3 and TM1 provide reasonable compression modulus 
(Kwr3 = 271.7 MeV, Krmı = 281.16 MeV) and asymmetry energy (Jura = 37.42 MeV, Jrmı = 36.89 MeV) 
but fairly small baryonic saturation density Lon = 0.145). In order to improve the description of these 
quantities, two nonlinear self-coupling effective interactions, PK1 with nonlinear o- and w-meson self-couplings 
and PKIR with nonlinear o-, w- and p-meson self-couplings were developed [51] (see Table 2). 

In order to reproduce better the experimental quantities such as the binding energies and nuclear radii, etc., 
an additional correction should be added in calculating the energy, i.e., the center-of-mass correction. Conven- 
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3 : 
tionally a phenomenological center-of-mass correction as = 41A~1/8 is used. Microscopically the correction 


can be calculated by the projection-after-variation in first-order approximation [87]: 


] 1 
Em = "OW (Pu (48) 


where the center-of-mass momentum P..m. = x p; and the expectation value of its square Pu) reads 


(Pim) = Kä vn, = > U2U5Pab g DA + > VaUaVbUbPab * Pab (49) 
a a,b a,b 


with occupation probabilities v2 and u2 = 1 — v2 accounting for the pairing effects, where a,b denote the BCS 
states (see the following section). 

Here it should be mentioned that the prescription (48) is based on non-relativistic considerations. It does 
not preserve Lorentz invariance. Furthermore, it also breaks the complete self-consistence of the variational 
scheme. However, as it is not included in the self-consistent procedure and only presents an additional correction 
term to the binding energy, one can be satisfied with it for the moment. Compared with the binding energy, 
this center-of-mass correction is sizable in light nuclei (about 9% in 160) but much less important in medium 
and heavy nuclei (about 0.4% in ?95Pb) as seen in Fig. 3. 


(MeV) 


E m. 


16 32 48 64 80 96 112 128 144 160 176 192 208 
Mass Number 


Figure 3: The microscopic center-of-mass correction in comparison with two phenomenological cases. Taken 


from Ref. [51]. 


For the nuclear radii, the effects from the center-of-mass motion can also be taken into account as follows. 
Because of its fairly small effects, a rather rough correction is adopted for protons, 


A Z A 
2 
ôR, = A. (al Re: 2 rila) +) (al Rem. Ida) » (50) 
where the center-of-mass coordinate Rem, = 1/A KE r;. Then one obtains, 
2 1 
2 2 2 
öR, = -gi + arm (51) 


where Rp and Rm denote the proton and matter radii. Here we only consider the direct-term contributions to 
Eq. (50) to conform to the spirit of the RMF theory. For the neutron radii, one can follow the same procedure 
as that for protons. The charge radius is obtained from the proton radius combining with the proton and 
neutron size and the center-of-mass correction (51) included in R2 [50] 


Hà, = RŽ + (0.862 fm)? — (0.336 fm)? N/Z. (52) 
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Table 2: The nonlinear effective interactions PK1, PK1R and density-dependent effective interaction PKDD 
[51] in comparison with TM1 [50], NL3 [80], TW99 [61], and DD-MEI [62]. 


PK1 PKIR PKDD TM1 NL3 TW99 DD-ME1 
Mn 939.5731 939.5731 939.5731 938 939 939 938.5000 
Mp 938.2796 938.2796 938.2796 938 939 939 938.5000 
Mo 514.0891 514.0873 555.5112 511.198 508.1941 550 549.5255 
Mu 784.254 784.2223 783 783 782.501 783 783.0000 
Mp 763 763 763 770 763 763 763.0000 
Io 10.3222 10.3219 10.7385 10.0289 10.2169 10.7285 10.4434 
Jw 13.0131 13.0134 13.1476 12.6139 12.8675 13.2902 12.8939 
Ip 4.5297 4.55 4.2998 4.6322 4.4744 3.661 3.8053 
g2 —8.1688 —8.1562 0 —7.2325 —10.4307 0 0.0000 
93 —9.9976 —10.1984 0 0.6183 —28.8851 0 0.0000 
C3 55.636 54.4459 0 71.3075 0 0 0.0000 
ds 0 350 0 0 0 0 0 


Table 3: Density-dependent parameters of PKDD [51] for meson-nucleon couplings in comparison with 
TW99 [61] and DD-ME1 [62]. 
Ag bs. Co do Ay) b, C d. ap 
PKDD 1.327423 0.435126 0.691666 0.694210 1.342170 0.371167 0.611397 0.738376 0.183305 
TW99 1.365469 0.226061 0.409704 0.901995 1.402488 0.172577 0.344293 0.983955 0.515000 
DD-ME1 1.3854 0.9781 1.5342 0.4661 1.3879 0.8525 1.3566 0.4957 0.5008 
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With the microscopic center-of-mass motion, a multi-parameter fitting can be performed using the Levenberg- 
Marquardt method [88]. In Ref. [51], the masses of 160, ^9 Ca, 45Ca, DON. 6®Ni, 907, 116Sn, 13295. !°*Pb and 
208Ph and the bulk quantities of nuclear matter are chosen as observables to determine the effective interactions. 
The radii are excluded because the proper values of the compression modulus K and the baryonic saturation 
density po are sufficient to give a good description of the radii. For a fixed value of the compression modulus 
K, a large baryonic saturation density po will give a small charge radius and vice versa. Therefore a proper 
description of both masses and radii of finite nuclei could be obtained by carefully adjusting the values of these 
two quantities K and po. To give a fairly precise description of the masses, the center-of-mass correction is also 
essential for both light and heavy nuclei. As it can be seen in Fig. 3, the deviation between the microscopic 
and phenomenological results is considerably large not only for the light nuclei but also for the heavy ones. 
And there exist very remarkable shell effects in the microscopic results which are impossible to obtain with the 
phenomenological methods. The microscopic center-of-mass correction [87], therefore, is chosen to deal with 
the center-of-mass motion. 

Because the contribution to the nuclear masses from the nonlinear p-meson term is found to be fairly 
small, the effective interaction PK1R is obtained by fixing the nonlinear self-coupling constant da to 350.0 and 
adjusting other parameters. 

In the RMF theory with density-dependent meson-nucleon couplings, the density-dependence of the coupling 
constants gg and g, can be parameterized as, 


Gi (Pu) = gi Pest) fiU) for i — 0,0, (53) 


V where 


an E, + b;(x KR di)? 
= 1+ Ci (a bi di)? 


>A is a function of r = py/Psat, and Psat denotes the baryonic saturation density of nuclear matter. For the p 


fi(z) (54) 


=) meson, an exponential dependence is utilized as 


Ip = gp Psat) exp[-ap(x — 1)]. (55) 


d „ For the functions Liz), one has five constraint conditions f;(1) = 1,f7(1) = f%(1) and f”(0) = 0. Then 


8 parameters related to density dependence for o-N and w-N couplings are reduced to 3 free parameters. In 


> ^ general, the masses of the nucleons and the p-meson are fixed and the nonlinear self-coupling constants g2, 93, c3 
(© and ds are set to zero. With 4 free parameters for density dependence, totally there are 8~9 parameters left 


free in the Lagrangian density (1) for the density-dependent meson-nucleon coupling RMF theory. A density- 
= dependent meson-nucleon coupling effective interaction PKDD has also been obtained in Ref. [51] (see Tables 2 
» and 3). 

Tables 2 and 3 tabulate the new effective interactions PK1, PKIR and PKDD [51] in comparison with the 
old ones TM1 [50], NL3 [80], TW99 [61] and DD-MEI [62]. The newly obtained ones reproduce better the 
experimental masses [89]. PK1, PK1R and PKDD also describe the charge radii very well, especially for those of 
the Pb isotopes. More comprehensive comparisons between Hartree-Fock-Bogoliubov, extended Thomas-Fermi 
model with Strutinski integral, RMF, and macroscopic-microscopic approaches with different forces have been 
performed for the description of nuclear masses and charge radii of spherical even-even nuclei (116 nuclides), 
from light (A=16) to heavy (A=220) ones in Ref. [90]. 

Table 4 lists the nuclear matter quantities calculated with the newly obtained effective interactions PK1, 
PKIR and PKDD, in comparison with those from the other interactions. All the new effective interactions 
give a proper value for the compression modulus K. 


2.4 Density and isospin dependence of effective interactions 


There are so far quite a number of effective interactions, PK1, PK1R, PKDD [51] together with NL1, NL2 [91], 
NL3 [80], NLSH [79], TM1, TM2 [50], GL-97 [92] and the density-dependent effective interactions TW-99 [61], 
DD-ME1 [62], etc.. It is very interesting to investigate the density and isospin dependence of the interaction 
strengthes of various effective interactions in the RMF theory and study their effects on nuclear matter [93]. 
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Table 4: Nuclear matter properties calculated with PK1, PK1R and PKDD [51] in comparison with those with 
TM1 [50], NL3 [80], TW99 [61], DD-MEI [62]. 


Interaction Psat. (fm?) Ey [MeV] K [MeV] J [MeV] M*/M (n) M*/M(p) 


PK1 0.1482 — 16.268 282.644 37.641 0.6055 0.6050 
PKIR 0.1482 — 16.274 283.674 37.831 0.6052 0.6046 
PKDD 0.1496 — 16.267 262.181 36.790 0.5712 0.5706 
NL3 0.1483 — 16.250 271.729 37.416 0.5950 0.5950 
TM1 0.1452 —16.263 281.161 36.892 0.6344 0.6344 
TW99 0.1530 —16.247 240.276 32.767 0.5549 0.5549 
DD-ME1 0.1520 —16.201 244.719 33.065 0.5780 0.5780 


Although for the nonlinear self-coupling effective interactions, the density dependencies are only embodied 
in the Klein-Gordon equations, it is still worthwhile to obtain a quantitative understanding of the coupling 
constants. 

In Fig. 4, the density dependencies of the interaction strengthes for g,(top), gu(middle) and g,(bottom) 
in nuclear matter as functions of the nuclear density are shown. The curves in the figures from top to bottom 
are labelled in the order of from left to right. The shadowed area corresponds to the empirical value of the 
saturation point in nuclear matter, i.e., Fermi momentum kp = 1.35 + 0.05 fm^! or density p = 0.166 + 0.018 
` fm~3. For the nonlinear effective interactions, the “equivalent” density dependence of the interaction strengthes 


> for c, w and p is extracted from the meson field equations according to [51]: 


Go ^ ge [9207 + 930°] /ps, (56a) 
Qs ~ Gu — Cau / pv, (56b) 
go ^ Ho deet VT / pa- (56c) 


ENI For the o-meson, TW-99 and DD-ME1 exhibit quite different behaviors for g5 compared with those of the other 


effective interactions in either magnitude or slope. In particular, the strengthes in TW-99 and DD-ME1 for 


> the density interval in Fig. 4 are almost two times larger than that of GL-97. Quite different results can also 
‘= be seen at the empirical nuclear matter densities. For the w-meson, except for TW-99, DD-MEI, TM1, and 


GL-97, the strengthes in all the other effective interactions are density-independent. However, the strengthes 
are closer to each other at the empirical saturation density than those of the o-meson although large differences 
can also be seen at low densities. For the p-meson which describes the isospin asymmetry, the strengthes for 
TW-99 and DD-ME1 show strong density dependence in contrast with those of the other effective interactions; 
while those of PK1, PK1R and PKDD are just in between [51]. 

In Fig. 5, the behavior of the binding energy per particle E/A as a function of the baryonic density p is 
shown for symmetric nuclear matter (left) and neutron matter (right). It is seen that all the density-dependent 
meson-nucleon coupling effective interactions give softer results than the nonlinear self-coupling ones , especially 
for neutron matter. The behaviors predicted by PK1, PK1R are much softer than that by NL3 and a little 
harder than that by TM1. The results from PKDD are slightly softer than that from DD-ME1 and much 
harder than that from TW99 at high densities. All these behaviors can be explained in the density-dependent 
meson-nucleon coupling framework. 

As what have been mentioned in expressions (56), the meson-nucleon coupling constants in the nonlinear self- 
coupling of mesons can be expressed as some kind of density-dependence. Fig. 4 shows the density-dependence 
of the coupling constants for the nonlinear self-coupling effective interactions and for the density-dependent 
version, where almost all the density-dependent coupling constants decrease with increasing density except for 
go of NL3, NLSH, NL1, which has strong o) self-couplings. On the other hand, the coupling constants go and gw 
of TM1, which has relatively weak c self-coupling (g3 = 0.6183) and strong w self-coupling (c3 = 71.3075), are 
smaller than the others, which means that TM1 provides relatively weaker scalar and vector potentials. This is 
the reason why TM1 presents the softer behavior than the other nonlinear self-coupling effective interactions. 
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Figure 4: Density dependence of the interaction strengths for ge (top), g,, (middle) and g,(bottom) in nuclear 
matter as functions of the nuclear density. The shadowed area corresponds to the empirical value of saturation 
point in nuclear matter (Fermi momentum kp = 1.35 + 0.05 fm^! or density p = 0.166 + 0.018 fm^?). Taken 
» from Ref. [51]. 


In Fig. 5, TW99 predicts the softest results because of its relatively small g, as compared with DD-ME1, 
PKDD and NL3, and large go as compared with PK1, PKIR and TM1 in Fig. 4. As it is known, the repulsive 
potential would be dominant at high densities. In Fig. 5, NL3 gives the hardest results because of its constant 
J and large g, even though its o-N coupling constant gg increases with the density. For PK1, PKIR and 
PKDD, which present the mid soft behaviors in Fig.5, the coupling constants also lie between the largest and 
the smallest in Fig. 4. There exists a significant difference between symmetric nuclear matter and neutron 
matter. The density-dependent effective interactions PKDD and DD-ME1 present similar trends as PK1, 
: PKIR and TM1 in symmetric nuclear matter but much softer in neutron matter, which may be interpreted 
by the density-dependence of g,. For the effective interaction PK1R, the density-dependence of the g, is fairly 
weak as compared with that of the density-dependent meson-nucleon coupling effective interactions. It can be 
explained by a very weak p-field, which generates neutron-proton asymmetry field. The behavior of g, with 
respect to the neutron-proton ratio is shown in Fig.6. As expected, the behavior is symmetric with respect to 
In(N/Z) and the density-dependence becomes more remarkable with the increase of the baryonic density and 
the neutron-proton asymmetry. 


3 Relativistic continuum Hartree-Bogoliubov theory 


Pairing correlations are due to the short range part of the nucleon-nucleon interaction and play important 
roles in open shell nuclei. À simple and commonly used method of dealing with pairing correlations is the BCS 
approach under the constant gap approximation. The BCS method can be combined easily with the relativistic 
mean field theory as shown in Refs. [21,78]. However, the conventional BCS method is not justified for exotic 
nuclei because it could not include properly the contribution of continuum states. Bogoliubov transformation 
is the generalization of the BCS scheme. By quantizing the meson field and making Gorkov factorization, the 
relativistic Hartree Bogoliubov (RHB) formalism was derived [95]. In order to describe the exotic nuclei, the 
relativistic continuum Hartree Bogoliubov (RCHB) theory was developed, in which the continuum states are 
discretized and the RHB equations are solved in the coordinate space [36, 42]. 

In this section, the pairing correlations, the conventional BCS approach, the continuum states, the resonant 
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Figure 5: The energy per particle E/A in symmetric nuclear matter (left) and neutron matter (right) as 
functions of the baryon density p. The shaded area indicates the empirical saturation region. The filled 
diamond presents the data taken from Ref. [94] as comparison. 
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Figure 6: The density-dependence of g, for PKIR with respect to the neutron-proton ratio N/Z. Taken from 
Ref. [51]. 
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BCS method and several methods for obtaining single particle resonant states are briefly reviewed. Finally, 
the detailed formalism of the RCHB theory together with the discussion on the effective pairing interactions 
are given. 


3.1 Pairing correlations 
3.1.1 BCS approximation 


For Fermion systems like atomic nuclei, Kramers degeneracy ensures the existence of pairs of degenerate, 
mutually time-reversal conjugate states, which could be coupled strongly by a short-range force, namely a 
pp-channel in the framework of mean field. Physically, such pairing correlations lead to the superfluidity. In 
Hartree approximation, the ground state properties could be described by filling the single particle levels from 
the bottom up to the Fermi level. The occupation probability of each single particle level is either zero or 
one. Due to pairing interaction, pairs of nucleons are scattered from the levels below the Fermi level to those 
above. Thus for the open shell nuclei, one has to deal with the occupation probabilities ranging from zero to 
one. Mathematically, this could be treated by introducing the concept of quasi-particles. For stable nuclei, the 
pairing gap can be extracted from the experimental odd-even mass difference and the BCS approximation is 
very simple and useful. For M-single particle orbitals in the Hartree approximation, for each particle orbital k 
(k = 1,--- , M), one can introduce the quasi-particle energy: 


Ey = Die — A + A2], (57) 


(C and the corresponding occupation probability: 


— variation: 


1 cr — A 
2 k 
— que 58 
#51. (58) 
M 
) where A is the Fermi level determined by the particle number N = > v? and A the pairing gap determined 
k=1 


T - from the experimental odd-even mass difference 
1 
A= 5 LEN + 2) - E(N +1) - (E(N + 1) - E(N))], (59) 


:— or deduced from certain pairing Hamiltonian Hpair when the experimental binding energies are unknown. For 
er monopole pairing interaction Hpair = G KR aj OF apap, there lies the pairing gap equation with classical 
4 kk’>0 


G A 
A=G) wr= vw (60) 
k>0 k»0 ® 


where u; satisfies u? = 1 — vz. For fixed values of the gap parameters A or pairing interaction strengthes G 
for neutrons and protons, the self-consistent solution is obtained by iteration. 

The BCS method can be easily implemented into the relativistic mean field description with the single 
particle orbitals in the relativistic Hartree level. However, it is only valid for bound states. For the extremely 
neutron-rich (or proton-rich) nuclei near the drip lines, the continuum states begin to contribute as schematically 
shown in Fig. 7. The conventional BCS approach will involve unphysical states and lead to the divergence. 
For these nuclei, one must either investigate the detailed properties of continuum states and include the 
coupling between the bound state and the continuum by extending the BCS method to resonant BCS theory 
(see Section 3.2.2) or generalize the BCS model by unifying the Hartree (or Hartree-Fock) scheme and the 
Bogoliubov transformation (see Section 3.3). 


3.2 Continuum states 
3.2.1 Single particle resonant states 
For stable open shell nuclei, the pairing gaps are of the order of MeV [16]. There are some evidences that 


the pairing gaps will increase towards the drip lines. For nuclei close to drip lines, the one or two particle 
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Figure 7: The schematic picture for the difference in the pairing correlation of stable and drip line nuclei. 
Normally the last filled nucleon in stable nuclei has 8 MeV binding, while it is near the threshold in drip line 
nuclei. The pairing correlation in exotic nuclei provides the possibility to scatter valence nucleons back and 
forth in the continuum. 


separation energies are usually very small. For example, the one proton separation energies 5, are 0.138 MeV, 
0.60 MeV and 0.29 MeV in ®B, !?N and #!Cl; the two neutron separation energies Son are 0.30 MeV and 1.33 
MeV in H Li and !^Be [89]. A small particle separation energy means that the Fermi level is very close to the 
threshold. Therefore, pairing correlations can scatter valence nucleons between bound and continuum states as 


— seen from Fig. 7. An unphysical continuum state distributes over the whole space while a physical resonance 
| S, state with positive energy stays inside the nucleus for a long time and has better asymptotic behavior. For 
f 3 the cases where the pairing correlations involve the continuum, such as exotic nuclei, it has been demonstrated 
D that the largest contribution to the pairing correlations is expected to come from the resonant continuum 
cO part [36, 37, 41, 42,96, 97]. Therefore it's very important to investigate single particle resonant states in exotic 
— nuclei. 


Various techniques have been developed to study resonant states in the continuum such as the R-matrix 


— theory [98] and the extended R-matrix theory [99], the K-matrix theory [100] and the conventional scattering 


j theory [101]. By discretizing the continuum, the contribution of the resonant states can be self-consistently 


== taken into account via a Bogoliubov transformation in coordinate space [36, 37, 42]. The bound-state-type 


.— methods have also been developed, including the real stabilization method [102], the complex scaling method 
ƏL (CSM) [103], and the analytic continuation in the coupling constant (ACCC) method [104]. Within the 


(U framework of the non-relativistic Schrödinger equation, the ACCC method has been applied to investigate 


the energies and widths for resonant states in light nuclei combined with few-body methods [105, 106], and 
to study single-particle resonant states in spherical and deformed nuclei by solving the Schródinger equation 
with Woods-Saxon potentials [107]. Recently the ACCC method and phase shift scattering method have been 
combined with the RMF model. 

The philosophy of the ACCC method is that a resonant state becomes a bound one if one increases the 
attractive potential. T'hen the energy, width, and wave function of the resonant state can be obtained by an 
analytic continuation carried out via a Padé approximant (PA) from the bound-state solutions. Compared with 
other bound-state-type methods, the ACCC approach is very effective and numerically quite simple because 
many methods available for bound-state problems can be used and the PA for analytic continuation can be 
implemented in easily. 

By solving the spherical relativistic Hartree equations in a meshed box of size Ro self-consistently, one 
can calculate the ground state properties of a nucleus. The vector potential V(r) and the scalar potential 
S(r), energies and wave functions for bound states are also obtained. By increasing the attractive potential 
(V(r) 4- S(r)) to A(V(r) + S(r)) , a resonant state will be lowered and becomes bound if the coupling constant 
A is large enough. Near the branch point Ao, defined by the scattering threshold k(Ao) = 0 [104], the wave 
number k(A) behaves as 

il t pe. (61) 
u(r = Ao); | — 0. 
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These properties suggest an analytic continuation of the wave number k in the complex A plane from the 
bound-state region into the resonance region by Padé approximant of the second kind (PAII) [104] 


Cota tenn? +... +err” 


k e kN] m4 ————————M———————— 62 
2 e) \ aoe dr? ES 162) 


where z = VÀ — Ao, and ene, CL, Oh, do,...,dw are the coefficients of PA. These coefficients can be de- 
termined by a set of reference points x; and k(x;) obtained from the Dirac equation with A; > Ao, i = 
12,.,L - N +1. With the complex wave number k(A = 1) = k, + iki, the resonance energy E and the width 


T 
T can be extracted from the relation e = E — 2, (E,T € R) and k? = e? — M?, i.e., 


(M? + k2 — k2)? + 4k2k? + (M? + k? — k?) 


— M, 
2 


T 24/ (M? + k2 — k?)? + AK2k? — 2(M? + k2 — k?). (63) 
Based on the Schrödinger and Dirac equations, the stability and convergence of the energies and widths for 
single-particle resonant states with square-well, harmonic-oscillator and Woods-Saxon potentials have been 
æ investigated and their dependence on the coupling constant interval and the order of the PA have been dis- 

. cussed [108, 109]. Combined with the RMF theory, the ACCC method was employed to study energies and 


NJ widths for resonant states in stable nuclei 160 and *Ca in Ref. [110] and to explore the single-particle resonant 
= states, including resonance parameters and wave functions, in exotic nuclei in Ref. [111]. 


The continuation in the coupling constant can be replaced by the continuation in k plane along the k(A) 


x5 trajectory determined by Eq. (62) to the point kg corresponding to the wave number for Gamow state, i.e., 
— kg = kU^NI(A = 1) [104]. Similarly, the wave function p(kr,r) for a resonant state can be obtained by an 
^ analytic continuation of the bound-state wave function y(k;,7) in the complex k plane. One can also prove 
— that the wave function o(k,r) is an analytic function of the wave number k in the inner region r < Ro where 


j the Jost function analyticity dominates [104]. Therefore, the technique, which has been adopted to find the 


H complex resonance energy, has been used to determine the resonance wave function p(kr,r). Firstly, the PA 
— are constructed to define the resonance wave function at any point r in the inner region (r « Ro) [104] 


NU j= Pr(k,r) _ ao(r) +aı(r)k + ag(r)k? +... + ag(r)k" (64) 
d Qn(k,r) 1+ bi(r)k + bo(r)k2 +...+0n(r)kN "` 

» where the coefficients a;(r) (i = 0,1,..., L) and b;(r) (j = 1,2,...,.N) are dependent on r. These coefficients 

can be determined by a set of reference points k; and y(k;,r) obtained from the Dirac equation with A; > 

Ao, (i= 1,2,...,L+N+1). The resonance wave function p(kr,r) = qU^Nl (Kg, r) (r < Ro) can be extrapolated 
in this way. 

'The Dirac equation can be rewritten as two decoupled Schródinger-like equations for the upper and the 

lower components respectively [32]. For neutrons, the Schródinger-like equation for the upper component reads 


dia? (r) kilki + 1) 
dr? r2 


in the outer region where V(r) ~ 0 and S(r) ~ 0. Its solution is the well known Riccati-Hankel function 


GË (r) + (E? — M?)GË (r) — 0, (65) 


hE =n Etz) Zei (66) 


where Alz) = zj&(z) is the usual regular solution, ie. Riccati-Bessel function, and A,(z) = zn,(z) is the 
irregular solution, i.e. Riccati-Neumann function. The outer wave function is matched to the inner wave 
function at r = rm < Ro 


gl gni.) = C(kR) [Jx(krtm) + D(kn)ôs(ERrm)] ; ES 
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where C(kg) is the coefficient for matching and D(kg) = tan 6,(kr) with ö.(kr) the phase shift. It has been 
found that 6,(kz) is almost a constant when rm is large enough. Given the upper component, the lower 
component F,(r) can be calculated from the relationship between the upper and the lower components and the 
resonance wave function is finally normalized according to the Zel’dovich procedures [104]. The wave function 
for the neutron resonant state v1gg/2 in 60Ca with different matching points is given in Fig. 8 where one finds 
the convergence of the wave function with respect to the matching point rm when rm changes from 14 fm to 
18 fm. 
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] Figure 8: Real parts of the upper and lower components of the radial wave function for the neutron v1gg9/2 
> state in 9? Ca calculated by the ACCC method with different matching points rm. Taken from Ref. [111]. 


The resonant states can be also investigated in the scattering phase shift method |41, 97, 112] where the 
". RMF equations are solved with the scattering-type boundary conditions. At large distances, where both the 
scalar and the vector potentials are zero, the RMF radial equations can be written in the form: 
dG 2 k(k+1) 
ES? F (a = EE 3 )G = 0, (68) 
1 dG sk 
- Iu I9. 
E+M dr r 
where a? = E? — M?. These equations are suited for fixing the scattering-type boundary conditions for the 
continuum spectrum. They are given by: 


(69) 


G = Car|cos(6)j;(ar) — sin(d)nj(ar)], (70) 
Coin . . 
=EFM [cos(0) j1 (ar) — sin(ö)nı_ı(ar)], (71) 


where jı and n, are the Bessel and Neumann functions and ô is the phase shift associated to the relativistic 
mean field. The constant C is fixed by the normalisation condition of the scattering wave functions and the 
phase shift 6 is calculated from the matching conditions. In the vicinity of an isolated resonance the derivative 
of the phase shift has a Breit-Wigner form, i.e. 

dó(E) = T/2 (72) 

dE (E; — E)? +T?/4’ 

from which one estimates the energy and the width of the resonance. In the vicinity of a resonance the radial 
wave functions of the scattering states have a large localisation inside the nucleus. Close to a resonance the 
energy dependence of both components of the Dirac wave functions can be factorized approximatively by a 
unique energy dependent function [113]. As in the non-relativistic case [114], this energy dependent factor is 
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Figure 9: Energies and widths for the neutron states v3p3/2, v3pi3, V2f7/2, V2fs/2, Vlhgj2, and vii in 
12277. Solid circles represent the results of the ACCC method, while open circles denote the results of the 
scattering method. Taken from Ref. [111]. 


^ the square root of the Breit-Wigner function written above, or, equivalently, the square root of the derivative 


CO of the phase shift. 
3 The energies and widths for the neutron resonant states v3p3/2, v3pi/5, V2f7/2, V2f5/2, V1hgj2, and v1%43/2 


= in 1??Zr are shown in Fig. 9. The results of the ACCC and the scattering methods are in good agreement with 
se) each other for most of these states. From both methods, v2f5/ and v1i;3/5 have large widths, while Gift 


^ and v1hg ja are very narrow. The width for v3pa;? from the scattering method is slightly larger than that from 
the ACCC calculation. For the resonant state v3p,/2, neither the energy nor the width can be extracted from 


— the scattering calculation. From a quantum mechanical point of view, these single-particle resonant states are 
A; quasi-stationary ones captured by centrifugal barriers. The decay width for a resonant state can be roughly 
‘= = explained by the penetration through the barrier. For those states with the same I, i.e., the same centrifugal 


barrier, the higher state has a larger width, e.g., for the two f states v2f5/2 and v2f7/2. Although v1hg/2 is 


> £ higher than v2 f5/2, its width is smaller because of higher centrifugal barrier. A similar argument also holds for 
“© a much higher but narrow state v1713/5. 


` 5 3.2.2 BCS approximation with resonant states 


The conventional BCS method is incapable of describing weakly bound nuclei due to the oscillating asymptotic 
behavior of single particle wave functions in the continuum [37]. Recently, many authors demonstrate that by 
picking up only those low-lying resonant states in the continuum, the BCS method can be successfully applied 
to weakly bound nuclei [41,97,115-117]. This is because the resonant states are well localized inside the nucleus 
and there is a large region outside the nucleus where the resonant wave functions have values close to zero 
before they start oscillating [97]. Therefore,the wave functions of the resonant states can be taken to be zero 
beyond a cutoff radius [97]. Meanwhile a realistic pairing potential, such as a zero-range 6-force, prefer to pick 
up the resonant states, whose wave functions have a large overlap with those of bound states below the Fermi 
surface, rather than the continuum [115,117]. In such a way, the HF+rBCS and RMF+rBCS results do not 
depend sensitively on the cutoff radius [97]. 

The BCS method can be easily extended to study the width effects of resonant states. The extended BCS 
equations for a general (finite range) pairing interaction including the contribution of the resonant continuum, 
referred as the resonant-BCS (rBCS) equations [41,97], are : 


A; = KS Vgjjujvj T > Van | gu(e)u,(e)v,(e)de ? (73) 
j y ` 
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Figure 10: The rms neutron radii (left panel) and the two-neutron separation energies (right panel) of even Zr 
isotopes as a function of the mass number A. Taken from Ref. [41]. 


A = 25 Verve, + 3 Kugener d gv (duy (dvu (dde , (74) 
v! v! 


H 


N= A. + D gt (e)v2(e)de . (75) 


= Here A; are the gaps for the bound states and A, are the averaged gaps for the resonant states. The quantity 


2j, + 1 dó 
gle) = SE is the total level density and 6, is the phase shift of the state with the angular momentum 


(l,jv). The factor o Le) can take into account the effect of the width and it is approximately a delta function for 
a very narrow resonance. The interaction matrix elements are calculated with the scattering wave functions at 
resonance energies and normalized inside a volume where the pairing interaction is active (see last paragraph). 

The BCS equations in the conventional RMF+BCS formalism can be easily substituted by these rBCS 
equations (73,74,75). To make this implementation, one can solve the RMF equations with the scattering- 
* type boundary conditions or any other methods like the ACCC method described in the last subsection and 
^ find the energies, widths and wave functions for the resonant states. This approximation scheme, called the 
RMF--rBCS method, has been applied to study neutron-rich Zr isotopes [41]. It was demonstrated that the 
sudden increase of neutron radii close to the neutron drip line depends on a few resonant states close to the 
continuum threshold. Including these resonant states, the RMF-BCS calculations give practically the same 
neutron radii and neutron separation energies as the RCHB calculations (see Fig.10). 


3.3 Relativistic continuum Hartree-Bogoliubov theory 
3.3.1  Bogoliubov transformation and continuum spectra 


The contribution of the resonant states can be self-consistently taken into account in the Hartree-Fock Bo- 
goliubov (HFB) theory solved in the coordinate space [37]. The advantage of the general HFB theory is that 
the variation method based on quasi-particle transformation unifies the self-consistent description of nuclear 
orbitals and the BCS pairing theory into a single variation theory. In Fig. 11, schematic pictures for different 
bases in the HF, HFB and the canonical basis have been given. For a given many-body system, one can choose 
any M-dimension orthogonal and complete basis to diagonalize the Hamiltonian. That is shown in the first 
column of Fig. 11. If only HF approach are considered, the Hamiltonian is diagonal in the HF basis. We get a 
M dimension basis, with a Fermi surface which divides the full open orbitals and full occupied orbitals. They 
are given in the second column. After the Bogoliubov transformation, we are working in the quasi-particle 
basis, which is enlarged into 2M-dimension. The quasi-particle orbitals are reflection symmetric with respect 
to the Fermi surface, as shown in the third column. They must be transformed into the physical space — 
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Figure 11: The schematic pictures of different basis. The first column is an arbitrary M-dimension orthogonal 
and complete basis to diagonalize the system Hamiltonian. The second column is the single particle levels in 
the HF space. The third column is quasi-particle levels in the quasi-particle HFB space which is enlarged into 
2M-dimension. The quasi-particle orbitals has a reflection symmetry with respect to the Fermi surface. The 
physical orbital in the canonical basis is given in the forth column. 


— canonical basis [35,42], in order to have a similar explanation as the HF space in the second column. For the 
5 Bogoliubov transformation involving the continuum, as the Fermi energy is negative for a bound nucleus, only 
N those resonant states in the canonical basis with proper asymptotic behavior could be picked up. In the forth 
> column, a schematic picture for the orbitals in the canonical basis is given, now the orbitals could be partially 
> „ occupied, which is another essential difference between the HF basis and the canonical basis. 


` 3.3.2 The formalism for the RCHB theory 


In the BCS approach, the gap equation takes into account only pairing between time-reversed continuum 
— States. A more general pairing between continuum states at neighboring energies can be taken care of by a 
© continuum Bogoliubov approach. The improvement to the BCS approximation is to introduce the concept 
EN of quasi-particle by the Bogoliubov transformation. Instead of introducing a mapping factor, the Bogoliubov 
ch transformation transforms the equation of motion into quasi-particle space. In the following we shall consider 
— the transformation from the single particle basis to the quasi-particle basis. The single particle operators in 
^. coordinate space ay,ay, are connected to the operators in the quasi-particle basis, Bi, Bl , i = LM, by the 
— following transformation: 


BN 3 s 3 Uir Vir (4) 
(ir) = [em aj = fe Vir Uir aj m 


D 


where i and r represents the row and column indices respectively. The basis of the quasi-particle states, Wy,;(r) 
and wy, (r), are defined via the single particle coordinate basis of ö-function as: 


e / 3 4 1 ( m) 

i = | dr ôfr-r es TT 

Ké e Ae > 
In spherical case, the Dirac spinor wave functions Yy;(r) and vy;(r) are similar to Eq.(20) and satisfy the 

following normalization relation: 


f ËrWW =T. (78) 


The variation method based on quasi-particle transformation can unify the self-consistent description of nu- 
clear orbitals, as given by the HF approach, and the BCS pairing theory into a single variation theory. As the 
equations of motion are self-consistently solved in the quasi-particle space, the convergence is guaranteed auto- 
matically. The unphysical continuum are excluded and the contribution from the resonance states with positive 
energies can be taken into account. From now on, the concept of continuum in Bogoliubov transformation is 
just the resonance states with positive energy as the unphysical continuum are excluded already. 
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Following the standard procedure of Bogoliubov transformation, a Dirac Hartree-Bogoliubov equation could 
be derived and the unified description of the mean field and pairing correlation in nuclei could be achieved. 
Using Green’s function techniques it has been shown in Refs. [21,95] how one can derive a relativistic Hartree- 
Fock-Bogoliubov theory: after a full quantization of the system the mesonic degrees of freedom are eliminated 
and, in full analogy to the non-relativistic case, the higher order Green’s functions are factorized in the sense 
of Gorkov [118]. Finally, neglecting retardation effects and the Fock term, as it is mostly done in relativistic 
mean field theory, one ends up with RHB equations as the following: 


x 20 DS Kai (79) 


h(r,r') = [a-p+V(r) + 8(M + S(r))] (80) 


is the Dirac Hamiltonian and the Fock term has been neglected. The pairing potential is 


where 


Agp (r,r) = =f r f PA Y Vav ge Eri riri riyen), (81) 


kk! 


where k, H. k, or k’ represent the other quantum numbers which together with the coordinate specify the single 


2 particle states. It is obtained from the pairing tensor, 


Kee (r, Y) = (lak iak sl) = Kä ver) ve" r^), (82) 


E;>0 


and one-meson exchange interaction Vis gi er’ ;rjrj) in the pp-channel. The nuclear density is given as 


Pe) p e). (83) 


k,E;>0 


In the RHB theory, the ground state |W) of the even particle system is defined as the vacuum with respect 


— to the quasi-particle: Oli = 0 for all i, i.e., |V) = [[; Bj|-), where |—) is the bare vacuum. For an odd 
„> particle system, the ground state can be correspondingly written as: |Y); = B IL j Bi|—), where j is the level 


‘© which is blocked. The exchange of the quasiparticle creation operator gi with the corresponding annihilation 
“= operator 8; means the replacement of the column j in the U and V matrices by the corresponding column in 


= the matrices V*, U* [35]. 


As the one-meson exchange interactions are not able to reproduce even in a semi-quantitative way the 
proper pairing in the realistic nuclear many-body problem, the interaction V in Eq. (81) is replaced by the 
phenomenological interaction which has been generally used in the pairing channel of the conventional HFB 
theory. As in Ref. [36] V used for the pairing potential (81) is either the density-dependent two-body force of 
zero range: 


1 T 
V(ri,r2) = Voö(rı = r2)] [1 = 0103] ( = e (84) 
P0 
with the interaction strength Vo and the nuclear matter density po or Gogny-type finite range force: 
V(ri,r) = A e(ri-r9/» (W, + BiP? — HiP” — MiP? P7) (85) 
i—1,2 


with the parameter ui, Wi, Bi, H; and M; (i = 1,2) as the finite range part of the Gogny force DIS [119]. 
The RHB equations can be solved in different bases as well. In many applications an expansion of the wave 
functions in an appropriate harmonic-oscillator basis of spherical or axial symmetry provides a satisfactory 
level of accuracy. For example, in Ref. [120] the RHB equations were solved by expanding the nucleon spinors 
U,(r) and V;(r) and the meson fields in a basis of spherical harmonic oscillators. However, for nuclei near the 
drip lines the expansion in the localized oscillator basis presents only a poor approximation for the continuum 
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states. Oscillator expansions produce densities which decrease too steeply in the asymptotic region at large 
distance from the center of the nucleus. The calculated rms radii cannot reproduce the experimental values, 
especially for halo nuclei. The transformed harmonic oscillator (THO) basis has also been employed in the 
solution of the RHB equations in configurational space [64]. Solving the RHB equations in Woods-Saxon basis 
is another alternative. 

In order to describe the coupling between bound and continuum states more exactly, the RHB equations 
and the equations for the meson fields should be solved in coordinate space [36, 42, 121] and one arrives at the 
RCHB theory. When spherical symmetry is imposed, the wave function can be conveniently written as 


GEC) GEO) 
y= "rd l i _ = | 
Mo AO uas | Vers ant. | me (86) 


The above equation (79) only depends on radial coordinates and can be derived as the following integro- 
differential equations: 


= o) Kur) - (E +A- V(r) + S¢)Fu(r) +r / Pane. ud 
un = ` Fy(r) +(E+A-—V(r) - S(r))Gu(r) +r / r'dr'A(r,r')Gy(r) = 0, , 
TE) V avr) + (E A Vt) = Se) ree) +r f rarat Eue) = 0, ui 
avin) - ER) - (G - X V(r) + S(r))Gv (r) e f Har AU d pep =. © 


) where the nucleon mass has been included in the scalar potential S(r). Instead of solving Eqs. (21) and (23) 


— self-consistently for the RMF case, now one has to solve Eqs. (87) and (23) self-consistently for the RCHB 


. case. The densities are calculated as 


(r) Mueve) - ëtt), 
Ant? py(r) lH), 
4mr^ps(r) = u mailt HF), 
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(88) 


Anr?pe(r 


e" = where the summations are over all quasi particle states. From the densities given above, one can calculate the 
(* rms radii and charge radius of the nucleus. The total binding energy is calculated as 


E = Enucleon + E, T Eu EE Ep + Ee ch Eco, (89) 


where except for Eynucleon, the other terms are the same as those in the spherical RH theory. The energy of 
nucleons Eyucleon is calculated as 


a = 2. 7 dr(A — E’) [|G (C)? + Ltr — 2Epair, (90) 


and Epair = -3TrAr. For the details of the transformation from the quasi particle space to the canonical 
basis, the reader is referred to Ref. [42]. 

To solve RCHB equations (87), one has to discretize the RCHB continuum. The discretization can be done 
by solving the RCHB equation in a spherical box of radius R, i.e., by imposing the boundary conditions 


Vu (£, |r| = R,s) = vy(E, |r| = R,s) = 0. (91) 


By increasing R, one can make the RCHB spectrum dense and better approximate the continuum, while 
increasing Emar allows one to take into account coupling to highly excited quasi-particle state. Both R and 
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Emax must be reasonably taken in such a way that the final results do not depend on the chosen R and Emaz- 
Only those solutions with E > 0 are necessary due to the special symmetry of RCHB equation [42]. If E is the 
solution of Eq. (79), —E must be another solution, but they are not independent. Both solutions are connected 
as following: 

vy(-E,r) = vy(E,r), vy(-E,r) = —dy(£,r). (92) 


So one can concentrate only on the solution with either positive energies or negative energies. 

Equations (87), in the case of ó-force Eq.(84), is reduced to ordinary coupled differential equations and 
can be solved with shooting method combined with Runge-Kutta algorithms, which is so far the most elegant 
method for coupled differential equations in coordinate representation. For each fixed value of the energy E 
the RCHB equations (87) are solved in the following steps: 


1. Discretizing the whole space [0, oo]; 
2. Choosing a value for the energy E in Eq. (87) with E > 0 for the given spin-parity channel; 


3. Finding the proper boundary condition for w(0), where ~ represents Gy,Fy,Gu,Fuy and integrated 
Eq. (87) outward from r = 0 to proper matching point r = Rmatch by Runge-Kutta algorithms; 


4. Finding the proper boundary condition for (oo) and integrated Eq. (87) from r = oo inward to proper 
matching point r = Rmaten by Runge-Kutta algorithms; 


5. Requiring the v(r) derived from the former two steps to be the same at r = Rmaten will lead to the 
correct energy E; 


6. Repeating the process till all the energies in the given spin-parity channel lying in |0, El have been 
found. 


© Equations (87) for Gogny forces (85) are a set of four coupled integro-differential equations. They are discretized 
^ in the space and solved by the finite element method as the following: 


1. Discretizing the whole space [0, oo] with N +1 points at r1 = 0,r2 = p,...,rN41 = N xp, where p = R/N. 


2. Assuming w(r) at rı = 0,r2 = p,..,TN41 = N x p are known, if p = R/N is small enough, v(r) in the 
element (i — 1) x p <r € i x p could be written as a simple linear combination of v(r;) and ıb(r;+1); 


3. Requiring Eqs. (87) is valid in the element (i — 1) x p € r € i x p, which could be easily integrated 
analytically and lead to a simple relation for E, v (ri) and w(ri+1); 


4. Repeating the process for all the N element, one then is lead to a 4 x (N + 1) linear algebra equations 
for E, Gy,Fy,Gp,Fpy for r = ri,T9,..., T N44, which is nothing but a generalized eigenvalue problem; 


5. Diagonalizing this generalized eigenvalue problem one can get all the energies and their wave function 
discretized at r = r1,T9, ..., T N41. 


Normally, equations (87) and (23) are solved in a self-consistent way by the shooting method or finite element 
methods with a step size of 0.1 fm using proper boundary conditions in a spherical box of radius R < 20 fm. 
The Fermi-surface A is determined self-consistently by the particle number [42]. The results do not depend on 
the box size for R > 15 fm [43]. For ö-force (84), the number of continuum levels is limited by a cut-off energy, 
which must be larger than the depth of the potentials. 


3.3.3 Pairing correlation: finite range vs. zero range 


In the vicinity of the drip line, as the contribution from the continuum is switched on, it is very challenging to 
have the proper single particle levels and treat the coupling of the continuum without any ambiguity. Therefore, 
a proper interaction for the coupling between the bound state and the continuum is also essential. In principle, 
the effective nucleon-nucleon interaction should be obtained by means of Bruckner renormalization which gives 
the correct interaction after modifying the free interaction for the effect of the nuclear medium. In practice, 
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Figure 12: Left: Two-neutron separation energies S2, of even Ni isotopes as a function of N, including the 
experimental data (solid points), RCHB with ó-force (open circles), RCHB with Gogny force (triangles) ‚HFB 
with SkP interaction (stars) and HFB with SIII interaction (pluses); Right: The neutron pairing energies for 
RCHB with Gogny force (triangles) and 6-force (open circles). The open circles connected by dashed line or 
solid line are respectively the pairing energies for RCHB with 6-force by fitting that of RCHB with Gogny force 
at ??Ni or Ni. Taken from [122]. 


C however, the effective interactions are approximated by some phenomenological forces, such as Skyrme-type 
J ö-force or finite range Gogny force. Delta force is relatively simple for numerical calculation and has a realistic 


density-dependent behavior, but it allows a coupling to the very highly excited states. Normally an energy 


— cutoff has to be introduced and the interaction strength has to be properly justified. The Gogny force has 
R * a better treatment for the coupling to the highly excited states, but it involves more sophisticated numerical 
= techniques. 


'The Skyrme- and Gogny-type pairing forces have been discussed before within the conventional HFB for- 


a malism, but the HFB equations for Gogny force are solved on a harmonic basis, which is valid only for stable 


nuclei. In Ref. [122], the proper form of the pairing interaction near the drip line is discussed in the framework 


^ of the RCHB theory. The chain of even-even nickel isotopes ranging from the proton drip line to the neutron 


- drip line are taken as examples. The pairing correlations are taken into account by either a density-dependent 


- force of zero range or the finite range Gogny force. Through the comparisons of the two neutron separation 
— energies Syn, the neutron, proton, and matter rms radii, excellent agreements have been found between the 


' calculations with both interactions and the empirical values. In Fig. 12, for example, the two neutron separa- 


tion energies Son of even Ni isotopes are shown as a function of the neutron number N from the proton drip 
line to the neutron drip line, including the experimental data (solid points), RCHB with ó-force (open circles), 
RCHB with Gogny force (triangles), HFB with SkP interaction [38] (stars) and SIII interaction [123] (pluses). 
The RCHB results with ó-force and Gogny force are almost identical. They show a strong kink at N=28 and 
a less one at N=50. The drip line nucleus is predicated at !"ONi in both calculations. The empirical data is 
known only up to N=50. Comparing with the available empirical data, the general trend and gradual decline 
of Sən have been well reproduced. 


The neutron pairing energies -3TrAk for the calculation of RCHB with Gogny force are given in Fig. 12 
as triangles. The pairing energies demonstrate the shell structure revealed by the So, in Fig. 12. They vanish 
at the closure shell and have a maximum values in the middle of two neighboring closure shells. A tendency of 
the enhancement of pairing energies appears for the neutron-rich nuclei. The pairing energies for RCHB with 
6-force are given by open circles. The circles connected by the solid line and dashed line are calculations with 
Vo = 742 MeVxfm~? (by fitting the pairing energy of the Gogny force at Ni) and 715 MeV xfm-^? (by fitting 
at ??Ni), respectively. The variation behavior of the pairing energies with N is quite similar, although Vo fitted 
at Ni gives 2 MeV more for stable nuclei and Vo fitted at ??Ni gives 4 MeV less for neutron rich nuclei. But 
these differences have little influence on the two-neutron separation energies and the rms radii, as one observes 
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Figure 13: The occupation probabilities in the canonical basis for HUN) as a function of the single particle energy 
around the threshold, for RCHB with ó-force (dashed line) and Gogny force (solid points). Taken from [122]. 


in Fig. 12, for example. This clearly demonstrated the reliability of the result against the pairing strength. 
Another interesting feature of the exotic nuclei is the contribution from the continuum. Definitely it is 

very interesting to see how different the contribution from the continuum is for the calculations with Gogny 

force and ó-force. In Fig. 13 the occupation probabilities of "Ni in the canonical basis is shown [35] for the 


“= neutron levels near the Fermi surface, i.e. in the interval —10 < E,,. < 10 MeV. The occupation probabilities 

“© for RCHB 6-force are represented by the dashed line and RCHB with Gogny force by the solid points. For 

> X Vo = 715 MeV xfm^? or Vo = 742 MeV xfm”?, the occupation probabilities are the same. Even the difference 

e < from the Gogny and the 6-force is too small to have any influence on the physical observables like Son and rms 
." radii. 


It has been shown that after proper renormalization (e.g., fixing the corresponding pairing energies), physical 


~. observables such as $2, and rms radii remain the same, even within a reasonable region of the interaction 


strength [122]. 


4 Exotic nuclei 


The study of exotic nuclei has attracted world wide attention due to their large N/Z ratios (isospin) and 


- interesting properties such as halo and skin. It is also important for other fields including astrophysics, e.g., 


— the properties of exotic nuclei are essential to understand the nucleosynthesis in the r process. Since the first 


case of halo in an exotic nucleus !! Li was observed with RIB in 1985 [13], more and more exotic nuclei have been 
investigated with various modern experimental methods to understand this attractive phenomenon better [6-8, 
11,12]. For nuclei far from the 6-stability valley and with small nucleon separation energy, the valence nucleons 
in exotic nuclei extend over quite a wide space to form low density nuclear matter. Furthermore, the Fermi 
surface for exotic nuclei is usually very close to the continuum threshold. The valence nucleons could be easily 
scattered to the continuum states due to the pairing correlation. T'hus, theories which can properly handle the 
pairing and continuum states are needed to describe the properties of exotic nuclei. One of such theories is the 
RCHB theory [36, 42] introduced in the previous section. In this section, the RCHB theory will be applied to 
study the properties of exotic nuclei, e.g., the binding energies, particle separation energies, the radii and cross 
sections, the single particle levels, shell structure, the restoration of the pseudo-spin symmetry, the halo and 
giant halo, and halos in hyper nuclei, etc. 


4.1 Binding and separation energies 


'The self-consistent microscopic RCHB theory has been used to study the ground state properties systematically 
such as binding energy, separation energy, radius, density distribution, single particle spectrum and shell 
structure etc. for a large number of isotope chains from the proton drip line to the neutron drip line throughout 
the periodic table including Li, C, N, O, F, Na, Ca, Ni, Zr, Sn and Pb [32, 36, 42-47, 122, 124-126]. In this 
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Figure 14: The two neutron separation energies $5, for the proton magic isotopes plotted as a function of 
the neutron number in RCHB with NLSH. The experimental values are denoted by the solid symbols and the 
calculated ones by the open symbols. Taken from Ref. [126]. 


subsection, we will briefly present the description of the binding energies and separation energies for some 
typical spherical nuclei in RCHB theory. 
Ground state properties of all the even-even O, Ca, Ni, Zr, Sn, Pb isotopes ranging from the proton drip line 


> to the neutron drip line were investigated with the RCHB theory [126]. The binding energy E, calculated from 
— the RCHB theory for these isotope chains and the corresponding data available [89] are compared. Generally 


speaking, the calculations reproduce the available data of binding energy quite well. The errors are usually 
' within 3-4 MeV which is less than 196 of the experimental values. Since the deformation effects are not included 
in the RCHB calculations yet, for some nuclei with large deformation, larger discrepancies (A Ey ~ 10 MeV) 
ə of binding energy from experimental values are found, e.g., in some Zr isotopes. 

In Table 5, the binding energies for Pb isotopes calculated with different effective interactions PK1, PK1R, 
. PKDD, TM1, NL3, TW99 and DD-ME1 are given in comparison with the experimental values [51]. Many 
effective interactions, particularly PK1, PKIR and PKDD provide good descriptions of the available data. All 


4 * the effective interactions overestimate the binding energy in the beginning of the isotopic chain. However, from 


. 190pp to ?10Pb. the effective interactions PK1, PKIR and PKDD give better descriptions than all the others. 
There exist however systematic deviations out of the neutron magic numbers, e.g., in ?'*Pb. For the effective 
interactions TM1, NL3, TW99 and DD-ME1, the deviations in ?!4Pb are smaller but fairly large in SPb. 

The one or two nucleon separation energies Sp(n) or 55,(5,) are quite sensitive quantities to test a microscopic 
theory. In Fig. 14 both the theoretical and the available experimental Son are presented as a function of neutron 
number N for the O, Ca, Ni, Zr, Sn, Pb isotope chains, respectively. The good coincidence between experiment 
and calculation is clearly seen. Seen along the S2, versus N curve in Fig. 14, strong kinks can be clearly seen 
at N = 8, N = 20,28,40, N = 28,50, N = 50,82 and N = 126 for O, Ca, Ni, Zr, Sn, Pb isotope chains, 
respectively. All these numbers correspond to the neutron magic numbers, which come from the large gaps of 
single particle energy levels. 

For the proton magic even mass nuclei studied in Ref. [126], the neutron drip line nuclei are predicted as 
280, Ca, Ni, 1367, 1785. and 7°Pb, respectively. In Refs. [127,128], the last bound nucleus in neutron 
rich Ca isotopes is predicted as "Ca with the HFB and Skyrme HF methods, respectively. The reason for that 
72Ca is bound is caused by the halo effect, which also results in the disappearance of the normal N = 50 magic 
number in the RCHB calculation as shown in Fig. 14 where no kink appears at N = 50 for Ca isotopes. In 
stable nuclei, the N = 50 magic number is given by the big energy gap between the 199,5 and 3s, 2d levels. 
However in the vicinity of the drip line region of Ca isotopes, the single particle energy of 3s1/2 would decrease 
with the halo effect and its zero centrifugal potential barrier, the large gap between 1gg/2 and 3s1/2 disappears. 
Therefore, the nucleus "Ca is no more a doubly magic nucleus and "?Ca is bound. This disappearance of 
the N = 50 magic number at the neutron drip line is due to the halo property. This is different from the 
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Table 5: Total binding energies of Pb isotopes (in MeV) calculated with the nonlinear self-coupling effective 

interactions PK1, PK1R and the density-dependent meson-nucleon coupling effective interaction PKDD [51], 

in comparison with the experimental values [89] and the results of TM1, NL3, TW99 and DD-ME1. 
A Exp. PK1 PKIR PKDD NL3 TM1 TW99 DD-ME1 
182 1411.650 1416.431 1416.619 1416.309 1415.184 1420.872 1417.202 1415.216 
184 1431.960 1435.548 1435.706 1435.478 1434.569 1439.768 1436.761 1434.569 
186 1451.700 1454.258 1454.382 1454.192 1453.511 1458.222 1455.879 1453.306 
188 1470.900 1472.586 1472.673 1472.504 1472.056 1476.272 1474.640 1471.639 
190 1489.720 1490.555 1490.603 1490.468 1490.247 1493.958 1493.109 1489.657 
192 1508.120 1508.197 1508.201 1508.116 1508.130 1511.317 1511.329 1507.411 
194 1525.930 1525.536 1525.494 1525.474 1525.733 1528.378 1529.309 1524.937 
196 1543.250 1542.592 1542.502 1542.545 1543.085 1545.163 1547.012 1542.262 
198 1560.070 1559.378 1559.236 1559.329 1560.199 1561.685 1564.338 1559.403 
200 1576.365 1575.893 1575.697 1575.831 1577.076 1577.942 1581.344 1576.365 
202 1592.202 1592.095 1591.843 1592.022 1593.679 1593.905 1598.084 1593.133 
204 1607.520 1607.851 1607.545 1607.770 1609.906 1609.477 1614.197 1609.676 
206 1622.340 1623.126 1622.765 1623.167 1625.725 1624.530 1630.122 1625.966 
208 1636.446 1637.443 1637.024 1637.387 1640.584 1638.777 1644.790 1641.415 
210 1645.568 1644.844 1644.345 1643.643 1647.969 1647.245 1650.888 1648.312 
212 1654.524 1652.155 1651.573 1649.873 1655.289 1655.549 1657.020 1655.174 
214 1663.298 1659.382 1658.718 1656.084 1662.551 1663.706 1663.196 1662.011 
[e 0.00033 0.00036 0.00042 0.00031 0.00068 0.00078 ` 0.00033 
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Figure 15: The one neutron separation energies Sn for the nuclei !17220,13724N,15726() and 17-25F by RCHB 
theory(open circles) and their experimental counterparts (solid circles). Taken from Ref. [46]. 


disappearance of the N — 20 magic number which is due to deformation. Recently a new magic number 
N = 16 has been discovered in neutron drip line light-nuclei region which is also considered as owing to the 


© halo effect [14]. It can be expected that near drip line, more traditional magic numbers would disappear and 


. new magic numbers would be found with the same mechanism. 

Systematic RCHB calculations with NLSH have been carried out for the C, N, O and F isotopes [46]. The 
_= one neutron separation energies Sp predicted by RCHB and the data [89] for the nuclei 7720, 137 N, 15-260 
~ and !77?F are shown in Fig. 15 as open and solid circles respectively. For carbon isotopes, the theoretical 
one neutron separation energies for 11715202207 are in agreement with the data. The calculated S, is less 
than 0 (-0.003 MeV) for the odd-A nucleus !?C which is bound in experiment. While for the experimentally 
^ unbound nucleus ?!C, the predicted value of S, is positive. From the neutron-deficient side to the neutron 
=a drip line, excellent agreement has been achieved for the nitrogen isotopes. As almost all relativistic mean field 
approaches do, the RCHB calculations overestimate binding for ?°?6O which are unstable in experiment. For 
fluorine isotopes, the Sn in '%?°-?9F are overestimated in contrast with the underestimated one in !$F. The 


‘© last bound neutron-rich nucleus is predicted as ??F. In general, the RCHB theory reproduces the S; data well 


considering that this is a microscopic and almost parameter free model. As discussed above, there are some 
discrepancies between the calculations and the empirical values for some of the studied isotopes. This may be 
due to deformation effects, which has been neglected in Ref. [46]. 


4.2 Radius, nuclear density distribution, and cross section 
4.2.1 Nuclear root mean square radii 


'The nuclear radii are important basic physical quantities to describe atomic nuclei. The calculated neutron, 
proton, matter, and charge radii for all even Ca isotopes against mass number A are plotted in Fig. 16. For 
nuclear charge radii re, the well-known parabolic behavior along 49749 Ca is not reproduced in RCHB calculation. 
It is mainly due to the improper spherical supposition to the real deformed 4244464 nuclei. It can be clearly 
seen that the radii rp, r4, Tm and re all increase with the mass number A. The mass radii rm as well as rn 
increase much faster than rp and re. In addition to the normal increase of r, and rm with A, a slightly up-bend 
tendency occurs in the neutron drip line region from A — 62 up to A — 72. 

For light nuclei, the RCHB calculations predict a slight increase of proton or charge radii when the proton 
drip line is approached [46]. The neutron and proton rms radii predicted by RCHB for the nuclei !1?7??C, 
12—24NŅ 14-260 and 16-25F are given in Fig. 17. The neutron radii for nuclei in each isotope chain increase 
steadily. While the corresponding proton radii remains almost constant for nuclei in each isotope chain except 
for proton rich ones. The increase of proton radii with neutron number decreasing in the proton rich side in 
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Figure 16: The rms radii for the proton, neutron, charge and matter distributions are plotted as a function 
of the neutron number for the Ca isotope. The available experimental values are denoted by solid symbols. 
Taken from Ref. [126]. 
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€ Figure 17: The neutron and proton rms radii predicted by RCHB for the nuclei 107?2C, 12724N.14726() and 


‘ !6-35F. Taken from Ref. [46]. 


; ` each isotope chain is responsible for the slightly increase in the charge changing cross section of these nuclei on 


' the target !?C as discussed in Section 4.2.4. 

In Fig. 18 the neutron radii r,, from the RCHB calculation for the even-even nuclei of the O, Ca, Ni, Zr, Sn 
and Pb isotope chains are ploted. The predicted r„ curve using the simple empirical equation r, = ro: NU? 
with ro = 1.139 normalizing to ?0®Pb is also represented in the figure. The simple formula for r, agrees with 
the calculated neutron radii except for two anomalies. One anomaly appears in the Ca chain above N = 40, 
and the other in Zr chain above N = 82. The increase of exotic Ni and Sn nuclei is not as rapid as that in Ca 
and Zr chains. The regions of the abnormal increases of the neutron radii are just the same as those for Son. 
Both behaviors are connected with the formation of giant halo (see Section 4.4 for details). 

It is interesting to investigate the neutron number dependence of the proton radii (or charge radii) and 
the neutron radii for these proton-magic isotope chains. In Fig. 19, the rms charge radii re obtained from the 
RCHB theory (open symbols) and the data available (solid symbols) for the even-even Ca, Ni, Zr, Sn and Pb 
isotopes are shown. The RCHB calculations reproduce the data very well (within 1.596). For a given isotopic 
chain, an approximate linear N dependence of the calculated rms charge radii r. is clearly seen. However, the 
variation of r. for a given isotopic chain deviates from the general accepted simple A!/? law (denoted by dashed 
lines in Fig. 19), which shows a strong isospin dependence of nuclear charge radii for nuclei with extreme N/Z 
ratio. Recently, these charge radii as well as the large amount of experimental data are investigated. A formula 
with the ZV3 dependence and isospin dependence is proposed which describes well the charge radii all over the 
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Figure 18: The rms neutron radii for the proton magic even mass nuclei plotted as a function of the neutron 
number. The solid curve corresponds to the radii calculated with the N!/3 formula. 


nuclear isotope chart, see subsection 4.2.3. [129]. 


4.2.2 Neutron skin thickness 


) The measurement of nucleon density distributions provides basic and important information in nuclear struc- 
` ture. The proton distribution of stable nuclei can be determined accurately by elastic electron or muon scatter- 
d ing. However, there is not comparable measurement of the neutron density distribution up to now. Although it 
= is difficult to obtain the neutron density distribution, the difference in radii of the neutron and proton density 
' distributions can be more easily determined [130]. Theoretically, the neutron skin thickness, rn — rp could place 
" important constraints on effective interactions used in nuclear models [131-133]. 
Many investigations have been done aiming at the description of the neutron skin thickness in the framework 
) of the RMF model. With the pairing correlation treated by a constant BCS approach, the neutron rms 


EN radii of some closed and open shell nuclei are investigated with the spherical RMF model [134]. It has been 


. demonstrated that with an increase in isospin of nuclei, the neutron rms radii in the relativistic mean-field 
approach with the parameter sets NL1 and NL2 are larger than both the empirical data and the predictions of 


e ` the Skyrme mean field [135]. This overestimation of the neutron skin thickness remains true for the effective 


interaction NL3, although the relativistic Hartree-Bogoliubov model with NL3+D1S give a good result of the 
neutron-skin thickness along the Sn isotopic chain compared to the available data [136, 137]. 

The overestimation of the neutron skin thickness is attributed to a poor description of the isovector channel 
of the nonlinear effective interaction in Ref. [62] where the RHB model is extended to include density-dependent 
meson-nucleon couplings. A density-dependent effective interaction DD-ME1 is proposed in which the improved 
isovector properties result in a better description of charge radii, and especially the calculated values of neutron 
thickness are in much better agreement with experimental data. As what is shown in Fig. 19, the difference 
between the values calculated with NL3 and DD-ME1 increases with the number of neutrons to about 0.1 
fm at N = 82, but then it remains practically constant for N > 82. While both interactions reproduce the 
isotopic trend of the experimental data, NL3 obviously overestimates the neutron skin. The values calculated 
with DD-ME1, on the other hand, are in excellent agreement with the experimental data. This result presents 
a strong indication that the isovector channel of the effective interaction DD-ME1 is correctly parametrized. 

In Ref. [133], the sensitivity of the neutron skin thickness S = ra — rp in ?5Pb to the addition of nucleon-o-p 
coupling corrections to a selection (PK1, NL3, S271, and Z271) of interactions in RMF model was investigated. 
With the higher order nucleon-o-p coupling, the PK1 and NL3 effective interactions lead to a minimum value 
of the neutron skin thickness S — 0.16 fm and the S271 and Z271 effective interactions yield even smaller values 
of S = 0.11 fm. The addition of these higher order correction terms to the conventional RMF models can be 
associated with some kind of multimeson exchange processes occurring in the inner (higher density) region of 
nuclei. In the inner region of a heavy nucleus such as SPb, one could imagine that multimeson exchange 
processes introduced through nucleon-o-p coupling corrections could play an important role in modifying the 
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Figure 19: DD-ME1 and NL3 predictions for the neutron skin thickness of Sn isotopes, compared with experi- 
mental data. Taken from Ref. [62]. 


short-range behavior of the nucleon-nucleon interaction, by softening the density dependence of the symmetry 
energy at high density (the inner part of nuclei) and thus leading to smaller values of neutron skin thickness S 
compared to previous effective interactions. 


) 4.2.8 Nuclear charge radii: isospin dependence and the Z!/* law 


^ Among all the size quantities describing nucleus, nuclear charge radius has been measured with very high 
=) precision by various techniques and methods experimentally. With its accuracy, the study of the nuclear 
charge radii is very important to understand not only the proton distribution inside the nucleus but also the 

exotic phenomena such as halo and skin observed in exotic nuclei. Particularly if one can get a simple and 
` reliable formula for nuclear charge radii, it will be very useful to extract the decoupling of proton and neutron 
: in the exotic nuclei and provide information for the effective nucleon-nucleon interaction widely used in all 
the nuclear models. The available charge radii data for A > 40 were examined and its global behavior was 
studied in Ref. [129]. Instead of the widely accepted A! law, a new Z!/? formula with isospin effect has been 


Co proposed. 


Based on the consideration of the nuclear saturation property, nuclear charge radius R. is usually described 
by the 4!/? law [16,35] 
Re = r44", (93) 


where A is the mass number and Re = 5/3 Te, with re the rms charge radius. For very light nuclei, because 
of their small A and large fluctuation in charge distribution due to the shell effect with short period, it seems 
that the charge distribution radius as a bulk property has little meaning. A detailed analysis of charge radii 
data for A > 40 shows that r4 is by no means a constant, but systematically decreases with A; i.e., ra ~ 1.31 
fm for light nuclei (A ~ 40) and ra = 1.20 fm for very heavy nuclei (see upper left panel in Fig. 20). This 
fact implies that some physics is missing in Eq. (93). 

Definite evidences of the violation of A!/? law have been found in the measurements of isotope shift in rms 
charge radii [138-140]. In particular, 6(r?) 442,4 values (associated with an addition of two neutrons) are often 
found to be considerably smaller compared to what is expected from the A!/? law (d(r?) 449.4 = 3x (r2) A) 
A typical example is that the observed charge radii of the calcium isotopes TI" "Ca remain almost the same 
(except a very little change induced by deformation or shell effect), though the mass number A has changed 
significantly. In contrast, there is also evidence that the observed ó(r?)4,5,4 values (associated with the 
addition of two protons) are often greater than what is expected from the AV? law (e.g., ó(r?) for 4Ti—“4Ca, 
50Ti—48Ca, etc.). 

Along the B-stability line, the ratio Z/A gradually decreases with A, i.e., for light nuclei Z/A ~ 1/2, and 
for the heaviest 8-stable nucleus 7225U, (Z/A)!/3 = 0.7285. It turns out that (1/2)!/3/0.7285 ~ 1.09, which is 


3T 


(a) Along the B—stability line (d) Along the B-stability line 


13 


R.=1.228A 


r/1.631 [fm] 


r,/1.228 [fm] 
5 


Figure 20: Comparison of descriptions of some formulas for nuclear charge radius Re. The parameters in the 
first four formulas are obtained through least square fitting to data on charge radii of nuclei with A > 40. 
Taken from Ref. [129]. 


very close to the ra ratio 1.30/1.20 shown in the upper left panel of Fig. 20. A naive point of view is that the 
, charge radius of a nucleus may be more directly related to its charge number Z, rather than its mass number 
` A. Therefore, compared to the A!/? law, a Z!/3 dependence for nuclear charge radii may be more reasonable, 


R, — rgZ1/3, (94) 


as noted in Refs. [138,139]. An analysis of the very limited data of charge radii then available showed that rz 
remains almost a constant, i.e., rz — 1.65(2) fm for A 7 40. Similar Z 1/3 dependent formula for nuclear charge 
radii was also proposed in Refs. |141, 142). With the physics for such a simple Z 1/3 Jaw open, one notes here 
) that the Z!/? law of Re includes an isospin dependence compared to the A!/? law as following: 


A N-Z\" r Ng 
Re rzZ TZ (5 7 ) A (1 A ) 


Tz uni 1N-Z 
aA (1 >): (95) 


A 


One can get ra = rz/2'/? = 1.31 fm and 0.33 for the coefficient of the isospin term. Equation (95) is 
approximately equivalent to that obtained in Ref. [140], R. = 1.25A1/3 (1 — 0.2212), by fitting the data. 
Obviously Eq. (94) has merits that only one parameter is needed. 

The Z!/? dependence of nuclear charge radii was also used to modify the Coulomb energy term in the 
semi-empirical nuclear mass formula [143], and it was found that the agreement between the calculated and 
experimental results was improved. Moreover, the AT! law for the nuclear giant (monopole, dipole and 
quadrupole) resonance energy (x 1/R) also could be improved, if the A-1/3 dependence is replaced by a Z- 1/3 
dependence [144]. 

In the past two decades, a vast amount of new experimental information on the electromagnetic structure 
of nuclear ground states of many nuclei has become available with high accuracy. The vast amount of improved 
experimental results have been investigated by using the A!/? and Z!/? dependence [129]. The measured 
charge radii for 536 nuclei with A > 40 are shown in Fig. 20. The dependence of charge radii on the quadrupole 
deformation 6 has been taken into account for the rare-earth deformed nuclei. 

In the upper left and right panels of Fig. 20, the charge radii for the most stable 159 nuclei with A > 40 
along the £-stability line have been analyzed by using the A!/? and Z!/? dependence. In the middle left and 
right panels, the same has been done for the data of 536 nuclei with A > 40. Two significant features can be 
observed: (A) On the one hand, the agreement between the data and the calculated results using the Z!/? 
dependence is much better than that using the Al/? law, i.e., while there exists a global regular decrease 
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Figure 21: Left: The rms charge radii versus the neutron number N for even-even Ca, Ni, Zr, Sn, Pb isotopes. 
The RCHB calculation with ó-force is represented by open symbols, while the corresponding data is denoted 
by solid symbols. The dashed lines represent the predictions by the AV? law with ra = 1.228 fm. Right: 
The experimental (solid symbols) and RCHB predicted (dashed lines) coefficient rz = R./Z!/? for the nuclear 
charge radii as a function of isospin quantity 7 = N/Z in even-even Ca, Ni, Zr, Sn and Pb isotopes. An 
asymptotic behavior is drawn as a solid line. Taken from Ref. [129]. 


> of r4 with A, rz nearly remains constant (rz = 1.631(11) fm). The relative rms deviations ø for the gi 
` dependence are much less than those for the 4!/? law. (B) On the other hand, though the rms deviation for 


— the Z!/? dependence is significantly reduced, an isospin induced scattering of the data in the middle panels 
c » in Fig. 20 can be also observed compared with that in the top panels. In Refs. [38, 140,145, 146], the isospin 
=) effects has been considered based on the A!/? law. As the Z!/? dependence can describe the nuclear charge 
^^ radii much better than the AV3 law, the Z 1/3 dependence is a more reasonable starting point for describing 
=~ the isospin dependence of nuclear charge radii. 


To make the isospin dependent Z!/? formula for nuclear charge radii be based on more strong foundation 
and be also valid for exotic nuclei, the charge radii of nuclei far from the (-stability line up to drip lines are 


— needed. However, such data is not well available, an alternative is to require that the new isospin dependent 
s= ZV/3 formula should assort with a reliable and microscopic nuclear model such as the fully self-consistent and 


= microscopic RCHB theory. 
For even-even Ca, Ni, Zr, Sn and Pb isotopes, the rms charge radii r. obtained from the RCHB theory (open 


© symbols) and the data available (solid symbols) for the even-even Ca, Ni, Zr, Sn and Pb isotopes are given 


in Fig. 21. The RCHB calculations reproduce the data very well. For a given isotopic chain, an approximate 
` linear N dependence of the calculated rms charge radii re is clearly seen in Fig. 21, which shows that the 
) variation of ro for a given isotopic chain deviates from both the simple Z!/? dependence and the simple A!/? 
law (denoted by dashed lines in Fig. 21). Therefore, a strong isospin dependence of nuclear charge radii is 
necessary for nuclei with extreme N/Z ratio. 

In Fig. 21, the data and RCHB predicted rz = R./Z 1/3 for the proton magic isotopes are also presented as 
a function of n = N/Z. It is clearly seen that the coefficient rz increases linearly with n (except some deviations 
due to deformation or shell effects) and the slopes are nearly the same for these isotopic chains. The linear 7 
(or isospin Tz — (N — Z)/2) dependence of rz for an isotopic chain may be understood as the effect of the 
first order perturbation correction of nuclear wave function due to an isospin Tz dependent interaction [147]. 
Based on the analysis of data in the middle and upper panels of Fig. 20 and RCHB prediction in Fig. 21, the 
following isospin dependent Z!/? formula for nuclear charge radii has been proposed [129]: 


Re=rz Z [1 +b- 7), n=N/Z, (96) 


where n* is 7 = N/Z for the nuclei along the 8-stability line which can be directly extracted from the nuclear 
mass formula [16], rz = 1.631(11) fm as obtained in upper right panel of Fig. 20, and b = 0.062(9) obtained 
from the least square fitting. 

The analysis of the available data using Eq. (96) with rz and b thus obtained is displayed in the lower right 
panel of Fig. 20. The data are reproduced better by Eq. (96) than by Eq. (94). Description of the data with 
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an isospin dependent AV/3 formula, Re = 1.228 AV/3 [1 + 0.029(n — n*)] is also shown in the lower left panel of 
Fig. 20 for comparison. The parameter for the isospin term is also obtained by least square fitting of the data. 
It can be seen that this formula describes the data less successfully than Eq. (96). 

In Refs. [140, 145], a modified A!/? formula with additional ~ (N — Z)/A and ~ 1/A terms made some 
successes on describing nuclear charge radius. In addition, this formua was extended with ~ A^? terms [38], 
which are considered to account for the isospin dependence of charge radii in light nuclei. Based on the fact 
that R,Z -1/3 is approximately a constant, where R, is the rms radius of proton density distribution, a formula 
for nuclear charge radius was proposed in Refs. [38,141,142], Re = ((r, 21/3)? + 0.64)V/?. Detailed analysis of 
these expressions show: (i) The Z!/? dependence is a better approximation to available R, data than the AU? 
one;(ii) The simple Z 1/3 dependent formula (94) gives nearly the same result as the formula suggested in Refs. 
[141,142]; (iii) The new isospin dependent Z!/? formula is better than the two-parameter isospin dependent 
AV3 formula but less better than the three-parameter one. A better agreement of the isospin dependence Z1/3 
formula of charge radii can be achieved by fitting both rz and b simultaneously. While in Eq. (96), the same 
rz is used as that in Eq. (94), which are derived only from the most stable nuclei. Thus the second term in 
Eq. (96) mainly corresponds to the isospin dependence of charge radii for the nuclei away from the 8 stability 
line. It is expected that the modified Z!/? formula (96) will become more useful with more and more data 
obtained for the nuclei far from the 6-stability line. 


4.2.4 Density distribution and cross section 


. The measurement of interaction cross sections plays an important role in exploring the properties of exotic 


^; nuclei because it can provide information of the nuclear density distribution and nuclear radii. Based on 
: 3 the measurement of interaction cross section with radioactive beams at relativistic energy, novel and entirely 
— unexpected features appear: e.g., the neutron halo and skin as the rapid increase in the measured interaction 
^f? cross-sections in the neutron-rich light nuclei. Systematic investigation of interaction cross sections for an 
=> isotope chain or an isotone chain can provide a good opportunity to study the density distributions over a wide 
^ range of isospin. From the mean field point of view, the nucleon density distributions in an atomic nucleus 
— determine the mean potential. Since the single particle levels and the shell structure are directly determined by 


; the mean fields, the study of the isospin dependence of shell structure and the mean potential, which become 


H highly diffuse near the particle drip line due to the large diffuseness in the density distribution, is also crucial to 


understand properties of unstable nuclei. Particularly the study of the surface diffuseness will provide another 


X mean to understand exotic nuclei because the surface diffuseness can be strongly related with the spin-orbital 
‘© splitting which is mainly determined by the diffuseness in the mean field potentials and could be measured 


experimentally. 


" The density distribution and shell structure Systematic RCHB calculations have been carried out for 
all the nuclei in Na isotopes with mass number A ranging from 17 to 45 with parameter set NLSH [45]. The 
calculated binding energies £p and the interaction cross sections with the Glauber Model are shown in Fig. 22. 
The calculated binding energies E are in good agreement with the empirical values [89]. The resonance states 
of "Na and !5Na (with a positive Fermi energy) are exactly reproduced. !?Na is bound but unstable against the 
proton emission, reproducing the experimental observation. The neutron drip line nucleus has been predicted 
to be Na. The difference between the calculations and the empirical values for the stable isotopes is from the 
deformation [45]. 

The density distributions are examined and the relation between the development of halo and shell effect 
within the RCHB theory is also studied in detail for Na isotopes [45]. For example, the density distributions 
for both the proton and the neutron in Na isotopes are given in Fig. 23. As seen in the upper part of Fig. 23, 
the change of the neutron density is as follows: for the nucleus with less neutrons (N), the density at the 
center is low and it spreads only to some smaller distance. With the increase of N, the density near the center 
increases due to the occupation of the Zei level, so does the development of neutron radius. In the lower part 
of Fig. 23, the proton density shows different behavior. T'he surface is more or less unchanged because of the 
Coulomb Barrier, but with the increase of N, the density of the center decreases due to the slight increasing at 
the tail. This is due to the symmetry energy. But as the density must be multiplied by a factor Avr? before 
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Figure 22: Upper: The interaction cross sections or of “Na isotopes on a carbon target at 950A MeV: the open 
circles are the result of RCHB and the available experimental data (A = 20 — 23,25 — 32) are given by solid 
circles with error bars. The dashed line is a simple extrapolation based on the RCHB calculation for ?5-?! Na. 
Lower: Binding energies for Na isotopes, the convention is the same as the upper part, but the RCHB result 
for particle unstable isotopes are indicated by triangle. Taken from Ref. [45]. 
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Figure 23: The neutron (upper) and proton (lower) density distributions in Na isotopes. The same figures but 
in logarithm scale are given as inserts to show the tail part of the density distribution. Taken from Ref. [45]. 
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Figure 24: The same as Fig. 23 but for matter density distribution. The upper part is given in logarithm scale 
(the radius is labeled at the top of the figure) and the radius ro at which p(ro) = 1074 for different isotopes is 


‘© given as inserts. Taken from Ref. [45]. 


^^ the integration in order to give the fixed Z, the big change in the center does not influence the outer part of 
^ the proton distribution very much. 


The matter densities for the even Na isotopes are given in Fig. 24 where the shortest tail in the total density 


^ occurs for ??—?Na, the most stable ones. For either the proton or the neutron rich sides, the tail density 
NI increases monotonically. The tail (r > 10 fm) of the proton rich nuclei is mainly due to the contribution of the 
— proton and that of the neutron drip line nuclei is mainly due to the contribution of the neutron. Compared 
— with the neutron-rich isotopes, the proton distribution for nucleus with less N has higher density at the center, 


e lower density in the middle (2.5 < r < 4.5fm), a larger tail in the outer (r > 4.5 fm) part. 


It is found that there is interesting connections between the matter distribution and the single particle level 
` distribution (see Fig. 25): after ??Na, as it is a sub-closure shell for the 1d; /2; the neutrons are filled in the 
> 251/72 and Laus, So the tail of the density for 27Na is two order of magnitude larger than Na at r = 10 fm, 
while the tail of the density from ?’Na to ?! Na is very close to each other. But as more neutrons are filled 
in, the added neutrons are filled in the next shell 1f7/2, 2p3/2 and 2p,/2. So again two order of magnitude's 
increase has been seen from ?! Na to ??Na, and then a gradual increase after ??Na. So it becomes clear that the 
rapid increase in the cross section is connected with the filling of neutrons in the next shell or sub shell. In the 
inserts in Fig. 24: the radius rg at which p(ro) = 1074 fm^? is given as a function of the mass number to see 
the relation between the shell effect and matter distribution more clearly. It is very interesting to see a slight 
decrease of rg from proton drip line to Na. The tail of ?°Na is the smallest. From ?9Na to ?! Na, one sees a 
almost constant ro. After a jump from ?! Na to ??Na, a rapid increasing tendency appears again. 


In Figs. 25 , the microscopic structure of the single particle energies in the canonical basis [35,42] is given. 
In the left panel of Fig. 25, the single particle levels in the canonical basis for the isotopes with an even neutron 
number are shown. Going from A — 19 to A — 45 we observe a big gap above the N — 8, N — 20 major shell , 
and N = 14 sub-shell. The N = 28 shell for stable nuclei fails to appear, as the 2p3/2 and 2p4/5 come so close 
to 1f7/2. When N > 20, the neutrons are filled to the levels in the continuum or weakly bound states in the 
order of 1/7/5, 2p3/2, 2P1/2, and 1/55. In the right part, the occupation probabilities in the canonical basis of 
all the neutron levels below E = 10 MeV have been given for ?Na to show how the levels are filled in nuclei 
near the drip line. 
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Figure 25: Left: Single particle energies for neutrons in the canonical basis as a function of the mass number. 
'The dashed line indicates the chemical potential. Right: The occupation probabilities in the canonical basis 
for Na. Taken from Ref. [45]. 


The surface diffuseness and the spin-orbital splitting The RCHB theory has also been used to sys- 
tematically study the tin isotopes from proton drip line to neutron drip line [124]. The spin-orbital splitting 


Eid Eia 
2l 4-1 

— versus A for the whole isotope chain is given in Fig. 26 for the neutron spin-orbital partners (1d3/2, 1d5/2),(197/2; 1992); 
= (lino, 1413/2), (1391/2; 19372), (Ubu, 1/72), and (1hg/2, 1511/5), and the proton spin-orbital partners (1d3/, 1ds/2), 
and (fu, 1/75). It is very interesting to see that the spin-orbital splitting for the neutron and proton is very 

close to each other, at least for (1d3/5, 1d5/5), and (1fs/2, 1 f7/2) cases. The splitting decreases monotonically 
from the proton drip line to the neutron drip line for all the partners. To see the underlying reason of this 
behavior, it is very helpful to examine the origin of the spin-orbital splitting in the Dirac equation. The equa- 

tion of motion for the nucleon moving in scalar and vector potentials could be de-coupled and reduced for the 
upper component and the lower component, respectively. If it is reduced in the lower component, it will be 
related with another interesting topic — the pseudo-spin symmetry discussed in next subsection. As for the 

: spin-orbital splitting, the Dirac equation can be reduced for the upper component as the following: 


Eis = (97) 


B 1 d2M-V+S)d 
E- por n CR m 
dr E+2M-V+S dr dr 
ais) 8 1 EUM =T Tu) eios 
r2 E+2M-V+S$Sr dr ? 
= — (E+2M-V+S)(E-V-S)G®(r), (98) 
where 
—l—1 j=1+41/2 
R= ndo (99) 
l j =l— 1/2. 
The spin-orbital splitting is due to the corresponding spin-orbital potential 
1 kd(2M — V + S) 
KEE 100 
E+2M-V+S$Sr dr (100) 


with some proper normalization factor. It is seen that the spin-orbital splitting in RMF is energy dependent 


and they depends on the derivative of the potential 2M — V + S as well as the particle distribution. Therefore 


d(2M MS 
one can introduce the so-called spin-orbital potential: Vj; = MEE The Vj, for H9Sn, 140Sn, and 
T 


170Sn, are given in the left panel of Fig. 27. The Lo, for both proton and neutron are almost the same, as the 
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Figure 26: The neutron spin-orbit splitting Ej, — mnm Ye versus the mass number A in tin 


2l 4-1 
isotopes for neutron (1d3/2, 1d5/2), (197/2; 199/2), (li11/2, 1413/2), (1p1/2; 19372). (15/2. 1/7/2), and (1hg/2, 1h11/3) 
orbital and proton (1d3/5, 1d5/5), and (15/5, 17/5) orbital, respectively. Taken from Ref. [124]. 


potential V — S is a big quantity (~ 700 MeV), the isospin dependence in the spin-orbital potential could be 


\ neglected. Therefore the proton and neutron Vj, are the same in the present model. That is the reason why the 

© spin-orbital splitting for the neutron and proton is very close to each other in Fig. 26. From "Sn to 170Sn, the 

CN amplitude of Vis decreases monotonically due to the surface diffuseness. Neutron potentials V — S for !°°Sn, 

T d 11091 12085. 13085. 14097 15091, 16085, and !Sn are given in the right panel of Fig. 27 and its inserted figure 

= gives the radii Ro at which V—S = 100 MeV as a function of the mass number. As seen in the above equations, 

© the spin-orbital splitting is related with the derivative of the potential V — S . The surface diffuseness happens 
“for both the vector and scalar potential, as well as for V — S and V + S. 


For the decline of the spin-orbital splitting, it comes from the diffuseness of the potential or the outwards 


— tendency of the potential. The diffuseness of the neutron potentials V + S are given in the left panel of 
EN Fig. 28 for Sn, 120Sn, Sn, 10g Sn, 16085 and 170Sn. It is seen that the depth of the potential 


decreases monotonically from the proton drip line to the neutron drip line and the surface of the potential 


„> moves outwards. The inserted figure gives the radii Ro at which V + S = —10 MeV as a function of the mass 
^ number. The proton potentials V + S for 11985,7985, 13085, 140Sn,150Sn, 160Sn, and 170Sn are given in the 
e right panel of Fig. 28. With the increase of the neutron number, the proton potentials are pushed towards 


== outside as well due to the proton-neutron interaction. 


Interaction cross sections For medium and high energy reactions, the interaction cross section Cint can be 
calculated with the Glauber model which is based on the individual nucleon-nucleon collisions in the overlap 
volume of the high energy colliding nuclei. In the framework of this model, the nucleus-nucleus reaction cross 
section can be obtained from nucleon-nucleon reaction cross section. Consider the collision of a projectile 
nucleus P on a target nucleus T' when the nuclei P and T' are situated at an impact parameter b relative to 
each other, the total reaction cross section can be written as 


g= 2m | bab (1—T(b)), (101) 


where T'(b) is the transparency function describing the probability that the projectile will pass through the 
target without interacting. T'(b) is given by 


T (b) = exp ET J peter. zp)dbpazr pr(br,zr)dbrdzr t(b — bp + by)! , (102) 
where oyy is the nucleon-nucleon reaction cross section suitably averaged over the interacting n-n, p-p, and 


n-p pairs [148,149]. The profile function t(b) for the NN scattering is taken as a delta function for zero-range 
nuclear interaction. 
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Figure 27: Left: The derivative of the neutron potentials V(r) — S(r) for H?Sn, 4°Sn, and "?Sn; Right: The 
neutron potentials V(r) — S(r) for °Sn, !1985,708n, 13085, H0gy 15055, 16055, and 170Sn. In order to examine 
the surface diffuseness more clearly, the radii Ro at which V(Ro) — S(Ro) = 100 MeV has been given as an 
inserted figure. Taken from Ref. [124]. 
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Figure 28: Left: The neutron potentials V(r) + S(r) for ! Sn, H0Sn,120Sn, 13055, 140g, 19085. 1605n, and 
170Sn. In order to examine the surface diffuseness more clearly, the radii Ro at which V(Ro) + S(Ro) = —10 
MeV has been given as an inserted figure; Right: The proton potentials V(r) + S(r) for "°Sn,!?°Sn, 130Sn, 
(ioun 19095. 16085. and !Sn. In order to examine the surface diffuseness more clearly, the radii Ry at which 
V(Ro) + S(Ro) = —10 MeV has been given as an inserted figure. Taken from Ref. [124]. 
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Figure 29: The proton density distributions predicted by RCHB for the nuclei 107220, 127274N, 147260 and 
16-25F in logarithm scale. Taken from Ref. [46]. 


To compare the cross section directly with experimental measured values, the densities pp p(r) of the target 
12C and the Na isotopes obtained from RCHB (see Fig. 23) were used in the Glauber model calculation. The 
. cross sections for reaction of Na isotopes at 9504 MeV on !?C have been compared with the experimental 


í ~) values [150] in the upper part of Fig. 22. The agreement between the calculated results and measured ones are 

fine. The cross section below 22Na changes only slightly with the neutron number, which means the proton 
— density has played an important role to remedy the contribution of less neutron. From ??Na to the neutron 

90 drip line, a gradual increase of the cross section has been observed. After 32Na, although no data exist yet, a 

© relatively fast increase has been predicted. In short, the measured cross section shows similar behavior as that 
œ of RCHB. 


N| Charge changing cross section Systematic investigation of interaction cross sections for an isotope chain 

. or an isotone chain can provide a good opportunity to study the density distributions over a wide range 
.— of isospin [45,150]. However the contribution from proton and neutron are coupled in the measurement of 
e - interaction cross section. To investigate possible differences in proton and neutron density distributions, a 
- combined analysis of the interaction cross section and other experiment on either proton or neutron alone are 
— necessary. The charge-changing cross section which is the cross section for all processes which result in a change 
“= ofthe atomic number for the projectile can provide good opportunity for this purpose. In Ref. [151], the total 


" charge-changing cross section o for the light stable and neutron-rich nuclei at relativistic energy on a carbon 
target were measured. The charge changing cross section o.. are calculated by using the Glauber Model and 
the proton density of projectiles and total density of target provided by RCHB theory [46]. 

The proton density distributions predicted by RCHB for the nuclei 197220,12724N,14726() and !°-?5F are 
given in Fig. 29 in logarithm scale. The change in the density distributions for each isotope chain in Fig. 29 
occurs only at the tail or in the center part. The density is multiplied by a factor Ar? before an integration in 
order to get proton number or radii, therefore the change of density in the central region does not matter very 
much. What is important is the density distribution in the tail part. Compared with the cases of neutron-rich 
isotopes, the proton distribution of nuclei with less N has higher density in the center, lower density in the 
middle part (2.5 < r « 4.5fm), a larger tail in the outer part (r > 4.5 fm) which gives rise to the increase of rp 
and Gece for the proton rich nuclei as will be seen in following figures. 

With the proton density from RCHB calculations for the projectile used for pp(bp, zp) in Eq. (102), the 
total charge-changing cross sections oce of the nuclei 197220, 12724N, 14-260 and 1625 P on a carbon target at 
relativistic energy are calculated and given in Fig. 30. The results of the Glauber model are represented by 
open circles and available data [151] by solid ones with error bars. The agreement between the calculation and 
the measurement is very good. 


'The charge changing cross sections change only slightly with the neutron number except for proton-rich 
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Figure 30: The total charge-changing cross sections Gec of the nuclei 197220,1274N.14726() and 16725 p on a 
carbon target at relativistic energy. The open circles are the result of RCHB and the available experimental 
data are given by solid circles with error bars. Taken from Ref. [46]. 


nuclei. This indicates that the proton density plays an important role in determining the charge-changing cross 
` sections Cec. A gradual increase of the cross section can be observed towards the neutron drip line. However, 


NJ the big error bars of the data cannot help to conclude anything here yet. It is shown clearly that the RCHB 
=~) theory, when combined with the Glauber model, can provide reliable description for not only interaction cross 
— section but also charge changing cross section. From comparison between Fig. 30 and Fig. 17, one finds similar 
9© trends of variations of proton radii and of charge changing cross sections for each isotope chain which implies 
— again the important role that proton plays in determining the charge-changing cross sections. 


— 4.3 Single particle levels and pseudospin symmetry 


ENS concept of pseudo-spin is that the single particle orbitals with j = I+1/2 and j = (J+-2)—1/2 lie very close in 
“== energy and can therefore be labelled as pseudo-spin doublets with quantum number ñ% = n—1, l = 1-1, and 5 = 
^^ s= 1/2. This concept is originally found in spherical nuclei [28,29] , but later proved to be a good approximation 


“© in deformed nuclei as well [152-154]. It is shown that pseudo-spin symmetry remains an important physical 


= concept even in the case of triaxiality [155]. The origin of pseudo-spin is proved to be connected with the 


special ratio in the strength of the spin-orbit and orbit-orbit interactions [156,157] and the unitary operator 
` performing a transformation from normal spin to pseudo-spin space have been discussed [156-158]. However, 
it is not explained why this special ratio is allowed in nuclei. 

The relation between the pseudo-spin symmetry and the RMF theory was first noted in Ref. [159], in which 
Bahri et al found that the RMF explains approximately the strengthes of spin-orbit and orbit-orbit interactions 
found by non-relativistic calculations. More details have been given in Refs. [153,160], in which it was suggested 
that the origin of pseudo-spin is related to the strength of the scalar and vector potentials. Ginocchio first 
revealed that pseudo-orbital angular momentum is nothing but the “orbital angular momentum" of the lower 
component of the Dirac wave function [30]. 

In the Dirac equation of the nucleon, when the scalar potential S(r) and the vector potential V (r) are equal 
in amplitudes but opposite in sign, i.e., S(r) + V(r) = 0, or more generally, d|S(r) + V (r)|/dr = 0, there is an 
exact pseudospin symmetry in single particle spectra [30-32]. Under these conditions, there should be some 
special relations between the four components of the Dirac wave functions of the pseudospin doublets which 
have been used to test the pseudospin symmetry in realistic nuclei [161—164]. 

'The pseudo spin symmetry in deformed nuclei has also been extensively studied in the framework of the 
RMF model. The conditions for an exact pseudo spin symmetry and the relationships for the lower components 
of the Dirac eigenfunctions were investigated in Refs. [165-168]. The relationships between the upper and lower 
components of the two states in the pseudo spin doublets have been studied thoroughly for realistic deformed 
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+ 


ij=141/2 versus the binding energy E = 


[Ej zu + LUES. 
os ee for Zr and Sn isotopes. From left to right, the pseudo-spin partners correspond 


2l+1 
to (1d3/2, 251/2), (1f5/2, 23/2), (197/2, 2d5j2) and (2d3/5, 381/2), respectively. Taken from Ref. [31]. 


relativistic eigenfunctions [164, 169]. 
| There are also many other investigations on the pseudo-spin symmetry in the framework of RMF the- 
© ory [170-177]. A thorough review is certainly beyond the purpose of the present paper. In the following, the 
: ` pseudo-spin symmetry in RCHB theory will be briefly reviewed. 

From the radial Dirac equation for spherical nuclei, similar as Eq. (98), one can derive the radial equation 
for the lower components: 


d? 1 dV+S)d, u K(1— wi 1 «dV +S), 4 
nv ue a ee 1ESWSSe a Se 
= -(E+2M-V+S)(E-V -S)F*(r), (103) 


>< where k(x — 1) = JI +1). It is clear that one can use either Eq. (103) or equivalently Eq. (98) to get the 
` eigenvalues E and the corresponding eigenfunctions. In Section 4.2.4, Eq. (98) has been used to discuss the 
spin-orbital splitting in connection with the corresponding spin-orbital potential. If the Eq. (103) is used instead 
1 x d(V + S) 
` E-V-Sr dr 
E for the same / will degenerate. This is the phenomenon of pseudo-spin symmetry observed in Refs. [28, 29]. 
As shown in Ref. [43], the particle levels for the bound states in canonical basis are the same as those by solving 
the Dirac equation with the scalar and vector potentials from RCHB. Therefore Eqs. (98) and (103) remain 
the same in canonical basis even after the pairing interaction has been taken into account. 

In Eq. (103), the term which splits the pseudo-spin partners is simply the PSOP. The hidden symmetry 
for the pseudo-spin approximation is revealed as d(V + S)/dr = 0, which is more general and includes V + 
S = 0 discussed in [30] as a special case. For exotic nuclei with highly diffuse potentials, d(V + S)/dr ~ 
0 may be a good approximation and then the pseudo-spin symmetry will be good. But generally, d(V + 
S)/dr = 0 is not always satisfied in the nuclei and the pseudo-spin symmetry is an approximation. However, if 

1 &d(V 4- S) vil —k) 
Voss ar 


the relative magnitude of the pseudo centrifugal barrier (PCB) 


: ` and the pseudo spin-orbital potential (PSOP) term, , is neglected, then the eigenvalues 


|< | |, the pseudo-spin approximation will be good. Thus, the comparison of 


d. 
; pone 3 aN and the PSOP can provide us with 


some information on the pseudo-spin symmetry. In Ref. [31], with the scalar and vector potentials derived 
from a self-consistent RCHB calculation, the pseudo-spin symmetry and its energy dependence have been 


discussed in realistic nuclei. Afterwards, the pseudo-spin splitting for exotic nuclei in Zr and Sn isotopes has 
been studied [32]. 
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Figure 32: The comparison of the effective pseudo centrifugal barrier (PCB) (E — V — p (dashed 
&d(V +S). 
dr 
» arbitrary scale for d3j2 (upper) and s;/5 (lower) in 1207r. The dashed lines are a 1d3/2 and 251/2, and the 

dot-dashed lines are for 2d3/2 and 3s;;5. The inserted boxes show the same quantities, but the ordinate is 
„ magnified and the abscissa is reduced to show the behaviors of the effective PCB and the effective PSOP near 
the nuclear surface. Taken from Ref. [31]. 


lines and dot-dashed lines) and the effective pseudo-spin orbital potential (PSOP) ` (solid line ) in 


E 
The redefined pseudo-spin orbital splitting AE = 


LE i41/2 + (+ I)E 4a 
2] +1 

J in Fig. 31. In both isotopes, a monotonous decreasing behavior with the decreasing binding energy is clearly 
seen. The pseudo-spin splitting for 3s1/2 and 2d3/2 is more than 10 times smaller than that of the Zero and 
1daj5. As far as the isospin dependence of the pseudo-spin orbital splitting is concerned, the splitting in Sn 
isotopes gives a monotonous decreasing behavior with the increasing isospin. Particularly for 25;/; and 1d3/2 
partners, the pseudo-spin splitting in 170Sn is only half of that in °°Sn. Just as one expects [31], the pseudo-spin 
symmetry in neutron-rich nuclei is better. In Zr isotopes, although the situation is more complicated (e.g., 
the effect of the deformation which is neglected here), the pattern is more or less the same, i.e., a monotonous 
decreasing behavior with the decreasing binding energy and a monotonous decreasing behavior with the isospin. 
From these studies, the pseudo-spin symmetry remains a good approximation for both stable and exotic nuclei. 
A better pseudo-spin symmetry can be expected for the orbital near the threshold, particularly for nuclei near 
the particle drip line. 


-E 


A = ZH versus the average single particle 
+ 


energy E = for the bound pseudo-spin partners in Sn and Zr isotopes are plotted 


To understand why the energy splitting of the pseudo-spin partner changes with different binding energies 
and why the pseudo-spin approximation is good in RMF, the PSOP and PCB should be examined carefully. 
Unfortunately, it is very hard to compare them clearly, as the PSOP has a singularity at E ~ (V + $). As one 


H 


— 
is only interested in the relative magnitude of the PCB and the PSOP, the effective PCB, (E-V —S jus 
" 


d(V +5 
and the effective PSOP, nays) are introduced for comparison. They correspond to the PCB and the 
r 


e 
PSOP multiplied by a common factor E — V — S respectively. 


The effective PSOP does not depend on the binding energy of the single particle level, but depends on the 
angular momentum and parity. On the other hand the effective PCB depends on the energy. Comparing these 
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Figure 33: The comparison of the effective pseudo centrifugal barrier (PCB) (E — V — S)——_,— (dashed 


dV+S 
lines) and the effective pseudo-spin orbital potential (PSOP) (solid line ) multiplied by the square 


of the wave function F of the lower components in arbitrary scale for das (upper) and s,/5 (lower) in 1207r. 
Taken from Ref. [31]. 


H ` two effective potentials one could see the energy dependence of the pseudo-spin symmetry. They are given in 
CNI Fig. 32 for 51/2 (lower) and d3/2 (upper) of 1207r in arbitrary scale. 
j The pseudo-spin approximation is much better for the less bound pseudo-spin partners, because the effective 


- ^ PCB is smaller for the more deeply bound states. This is in agreement with the results shown in Fig. 31. The 
JO effective PSOP and the effective PCB are also given as inserts in Fig. 32 in order to show their behavior near 


— the nuclear surface. 


The effective PCB (dashed lines or dot-dashed lines) and the effective PSOP (solid lines) multiplied by the 
^, squares of the lower component wave function F (r) are given in Fig. 33, for 2s1/2 (upper left), 35; / (lower left), 
m ids /2 (upper right), and 2d3/, (lower right) of 1207r in arbitrary scale. The pseudo-spin approximation is much 


Z" better for the less bound pseudo-spin partners, because the effective PCB is smaller for the more deeply bound 


-== states. This is in agreement with the results shown above. It is clear that the contribution of the effective PCB 
^^ (dashed lines or dot-dashed lines) is much bigger than that of the effective PSOP (solid lines). Generally the 
i ` effective PSOP is two orders of magnitude smaller than the effective PCB. The energy difference between the 
= orbital j =i+1 /2 and the orbital j = TEN /2 is always negative with some exceptions. Further investigations 


£ À i d(V +5 ; er , 
» reveal that the integration of AV ie |? over r gives the splitting of the pseudo-spin partners, whose sign 


will decide the normal splitting or the reverse. 


4.4 Halos and giant halos 


Since the experimental discovery of neutron halo phenomena in ! Li [13], the study of exotic nuclei has become 
a very challenging topic in nuclear physics. Experimentally more and more halo nuclei have been observed [6- 
8,11,12,178-181]. On the theoretical side very different models have been used to describe these phenomena. 
Three-body calculations based on an inert core of ?Li and two outside neutrons treat the translational invariance 
properly, but they neglect polarization effects. Full shell model calculations contain all the configuration mixing, 
but they are limited to relatively small configuration spaces in the oscillator model. In the mean field description, 
the large radius observed experimentally in !! Li has been interpreted by the fact that filling in more and more 
neutrons in the nuclear well, the Fermi surface for the neutrons comes close to the continuum limit. Due to 
the small single-neutron separation energy, the tails of the two wave functions of the last filled orbital (1pı/2) 
reach very far outside of the nuclear well and a neutron halo is formed |13]. Although this simple interpretation 
is based on the mean field picture, several microscopic theoretical investigations within self-consistent mean 
field models [182-185] have failed. In fact, it would be a strange accident if the last occupied neutron level 
in HLi, the 1p, ya Orbit, would be so close to the continuum limit that the tail of the last two uncorrelated 
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neutron wave functions reaches so far out as would be necessary to reproduce the large experimental radius 
observed. Therefore, in the non-relativistic scheme Bertsch et al [182] and Sagawa [184] introduced an artificial 
modification to the potential in order to reproduce the small separation energy in their mean field calculations 
using Skyrme interactions. In this way the authors were able to reproduce qualitatively the observed trends, 
although some discrepancies still remained. 

Koepf et al [183] were the first to investigate the neutron halo in the RMF theory for both spherical and 
axially deformed cases. They found large deformations for the lighter Li isotopes, but a spherical shape for !! Li 
and as in the non-relativistic investigations the binding energy of the 1p;/2 was too large so as to reproduce a 
neutron halo. Zhu et al [185] improved the result of the RMF calculations by applying a similar modification 
to the potential as in Ref. [184] in order to adjust the proper size of the halo. 

The first self-consistent microscopic description of the halo in !! Li was achieved by using the RCHB theory 
where the pairing correlations are taken into account by a density dependent force of zero range and the 
halo was reproduced in a self-consistent way by using the scattering of Cooper pairs to the 2s;/2 level in the 
continuum [36]. 

In Ref. [121], pairing correlations and the coupling to particle continuum states are described by finite range 
two-body forces and the finite element method are used in the coordinate space discretization of the coupled 
system of Dirac-Hartree-Bogoliubov integro-differential eigenvalue equations, and Klein-Gordon equations for 
the meson fields |71]. The authors of Ref. [121] performed calculations for the isotopic chains of Ne and C 
nuclei and found evidence for the occurrence of neutron halo in heavier Ne isotopes. The RHB theory is also 
applied to the description of the properties of light nuclei in isotopic chains of C, N, O, F, Ne, Na, and Mg 


r . with large neutron excess in Ref. [186]. The RHB model predicts the formation of neutron skin and eventually 
( Y a neutron halo in Ne and Na due to the quasi-degeneracy of the triplet of states 17/5, 2p3/2, and 2p,/5 which 
=) are close to the Fermi surface in neutron rich nuclei. 


To describe medium-heavy nuclei far from the 6 stability line, the main problem of nuclear structure 


0 models is the extrapolation of effective forces to nuclei with extreme isospin values because most of these nuclei 
— are still not accessible experimentally and therefore one could not compare theoretical results with empirical 


data. In order to make predictions for medium-heavy nuclei at the neutron drip line, one has to test available 


— effective interactions in detailed calculations of the properties of neutron-rich nuclei for which a comparison 


; with experimental data is possible. This has been done in Ref. [187] where the RHB theory is applied in the 


H description of the ground-state properties of Ni and Sn isotopes. In a comparison with available experimental 


. data, it is shown that the NL3 + Gogny D1S effective interaction provides an excellent description of binding 


Ee energies, neutron separation energies, and proton and neutron rms radii. The results indicate that this choice 
* of model parameters might also be valid for nuclei with more extreme isospin values, i.e., medium-heavy nuclei 
— at the drip lines. 


We note that the RHB model solved in a discretized coordinate space is available only for spherical systems 
/ up to now. The deformed RHB theory solved in a harmonic oscillator basis [188] is not applicable for halo 
nuclei due to the localization property of the harmonic oscillator potential. One should work in two or three 
dimensional spatial lattice or instead, one can improve the asymptotic behavior of the HO wave function [64,65] 
or adopt other basis which has a correct asymptotic behavior [73]. 

In the following, we will briefly review the RCHB description of halo phenomena and the prediction of giant 
halo in light and medium-heavy nuclei. 


4.4.1 Neutron halo in light nuclei 


The ground state properties of Li isotopes with mass numbers A = 6 to A = 11 were investigated by using the 
RCHB theory with the effective interaction NL2 [36]. Pairing is neglected for the three protons and the strength 
Vo of the pairing force (84) for the neutrons is determined by a calculation in the nucleus "Li adjusting the 


corresponding pairing energy -~Tr Ax to that of a RHB-calculation in an oscillator basis using the finite range 


part of the Gogny force DIS of Ref. [119] in the pairing channel. Good agreement with experimental values is 
found for the total binding energies and the radii of the isotope chain 9Li to !! Li. The matter radius shows a 
considerable increase when going from the nucleus HL to !! Li. In contrast to the earlier mean field calculations 
of Refs. [182, 184,185] these results are obtained without any artificial modifications of the potential. 
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Figure 34: Calculated and experimental density distribution in !!Li and ?Li. The solid line shows the result 
of ''Li while the dashed line corresponds to the calculation of ?Li. The shaded area gives the experimental 
results with error bars. Taken from Ref. [36]. 
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Figure 35: The mean field potential S + V for protons (l.h.s.) and neutrons (r.h.s.). The chemical potential is 
- given by a dashed line. The energy levels in the canonical basis are indicated by horizontal lines with various 
lengths proportional to the occupation of the corresponding orbit. Taken from Ref. [36]. 


In Fig. 34 the corresponding density distribution for the neutrons in the isotopes ?Li and !!Li is shown. 
It is clearly seen that the increase of the matter radius is caused by a large neutron halo in the nucleus "Li. 
Its density distribution is in very good agreement with the experimental density of this isotope shown with its 
error bars by the shaded area. 

In order to understand the microscopic structure of this halo, in Fig. 35 the mean field S(r) 4- V(r) for the 
protons and neutrons together with the energy levels in the canonical basis is shown. The Fermi level for the 
neutrons is very close to the continuum limit in close vicinity to the v1p;/5 and to the v25;/5 level. The length 
of these energy levels in Fig. 35 is proportional to the corresponding occupation. One clearly finds that pairing 
correlations cause a partial occupation of both the v1p;/5 and the 725,5 level, i.e. a scattering of Cooper pairs 
to the continuum. 

In contrast to many previous investigations, in the RCHB model, the halo is not formed by two neutrons 
occupying the 1p1/2 level very close to the continuum limit, but is formed by Cooper-pairs scattered mainly 
in the two levels 1p,/5 and 2s;/5. This is made possible by the fact that the Zauuz comes down close to the 
Fermi level in this nucleus and by the density dependent pairing interaction coupling the levels below the Fermi 
surface to the continuum. The RCHB theory provides a much more general mechanism for halo occurring in 
nuclei. One needs only several single particle levels with small orbital angular momenta and correspondingly 
small centrifugal barriers close to, but not necessarily directly at, the continuum limit. 
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Figure 36: Canonical basis single-particle neutron levels as functions of the number of neutrons. The spectrum 
is calculated for Ne isotopes. Taken from Ref. [121]. 


In Ref. [121], the isotopic chains of Ne and C nuclei were studied with the RHB model with finite range 
pairing interaction in the coordinate space. The formation of neutron skin and eventually a neutron halo in Ne 
was found to be mainly due to the quasi-degeneracy of the triplet of states 1f7/2, 2p3/2, and 2p1/2. In Fig. 36, 
the neutron single-particle spectrum in the canonical basis is displayed for Ne isotopes. The energies of levels 


© in the continuum decrease with increasing neutron number (cf. Fig. 41 (a)). The shell structure dramatically 


` changes at N > 22. The triplet of states 1f7/2, 2p3/2, and 2p,/2 approaches zero energy, and a gap is formed 
| between these states and all other states in the continuum. The Fermi level uniformly increases toward zero for 
= N = 22. Between N = 22 and N = 32, the Fermi level is practically constant and very close to the continuum. 


^" The addition of neutrons in this region of the drip does not increase the binding. Only the spatial extension of 


. neutron distribution displays an increase. At N = 32 the Fermi energy becomes slightly positive, and heavier 
isotopes are not bound any more. The finite range pairing interaction promotes neutrons from the 1f7/2 orbital 
> to the 2p levels. Since these levels are so close in energy, the total binding energy does not change significantly. 


C ` Because of their small centrifugal barrier, the 2p3/; and 2p,/2 orbitals form the halo. 


> < 4.4.2 Giant halos in Zr isotopes 


In all of the halos observed so far, one has only a very small number of nucleons, namely one or two outside 
of the normal core. In order to study the influence of correlations and many-body effects it would be very 
interesting to find also nuclei with a larger number of neutron distributed in the halo. Following the successful 
description of the halo nucleus ‘Li, the RCHB theory was used to study and predict giant neutron halos for 
Zr-isotopes close to the neutron drip line [43]. It is formed by up to six neutrons outside of the !??Zr core with 
the magic neutron number N = 82. 

From the upper panel of Fig. 37 which presents the rms radii of the protons and neutrons for the Zirconium 
isotopes, one finds a kink for the neutron rms radius at the magic neutron number N = 82. This kink can 
be understood more clearly by considering the microscopic structure of the underlying wave functions and the 
single particle energies in the canonical basis which is shown in the lower panel of Fig. 37. Going from N = 70 
to N = 100, there is a big gap above the (hu orbit. For N > 82, the neutrons are filled to the levels in the 
continuum or weakly bound states in the order of 3p3/2, 2f7/2, 3P1/2, 2fs5/2 and Lhg/z. 

The neutron chemical potential is given in this Figure by a dashed line. It approaches rapidly the continuum 
already shortly after the magic neutron number N = 82 and it crosses the continuum at N = 100 for the nucleus 
to Zr. In this region the chemical potential is very small and almost parallel to the continuum limit. This 
means that the additional neutrons are added with a very small, nearly vanishing binding energy at the edge of 
the continuum. The total binding energies E for the isotopes above !??Zr are therefore almost identical. This 
has been recognized already in Ref. [189] in RMF calculations using the BCS approximation and an expansion 
in an oscillator basis, which is definitely not reliable for chemical potentials so close to the continuum limit. 

In the upper panel of Fig. 38 the density distributions for neutrons and protons in the nucleus Zr are 


53 


RMS Radius 


s.p.energy[MeV] 


Figure 37: Upper: rms radii for neutrons and protons in Zr isotopes close to the neutron drip line as a function 
of the mass number A. Lower: single particle energies for neutrons in the canonical basis. The dashed line 


indicates the chemical potential. Taken from Ref. [43]. 
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Figure 38: Upper: neutron and proton density distribution in !?^Zr. Dashed lines indicate calculations for 
different values of the box size R and the dashed-dotted line give the neutron distribution for the core (Zär. 
Lower: relative contributions of the different orbits to the full neutron density as a function of the radius. The 
shaded area indicates the total neutron density in arbitrary units. Taken from Ref. [43]. 
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` reliably reproduced only up to r = 25 fm, where it has decreased to 10-9 fm 
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Figure 39: The occupation probabilities in the canonical basis for various Zr isotopes with mass number A as 
a function of the single particle energy. The chemical potential is indicated by a vertical line. For A > 122 the 
number N, of neutrons in the halo is also shown. Taken from Ref. [43]. 


) shown. A very large box size is needed to describe the halo properly. For R —30 fm the neutron density is 
"2. The full line in the upper 
panel of Fig. 38 is an asymptotic extension to infinite box size. On the other hand, the density distribution 
inside the nucleus is reproduced properly even for small values of the box size. The relative contributions pj; 
of the different orbits characterized by the quantum numbers nlj with respect to the total neutron density p; 
„ are displayed in the lower panel of Fig. 38. The shaded area gives the total neutron density in arbitrary units. 
The halo is formed essentially by contributions from three orbits 3p3/2, 3p1/2, and 2f7/2. The most inner part 
of the halo (7 < r < 9 fm) the 2f7/2 orbit plays the dominant role. As can be seen in the lower part of Fig. 38, 
this orbit is slightly below the chemical potential and to the continuum limit in this nucleus. Further outside 
. (10 € r € 15 fm) its relative contribution is strongly reduced because of the larger centrifugal barrier felt by 
' the 1 = 3 orbit. In this region the orbit 3p; /2, Which has nearly the same position as the 2f7/2 orbit, takes 
over. Because of the smaller orbital angular momentum, this orbit feels a reduced centrifugal barrier. For even 
larger distances from the center (r > 15 fm) its relative contribution is somewhat reduced and the 3p; j? orbit 
gains importance. The 3p3/2 and the 3pj/2 levels feel the same centrifugal barrier, but the latter is situated 
directly at the continuum limit and therefore it is more loosely bound than the other two orbits. 


'The occupation probabilities in the canonical basis of all the neutron levels near the Fermi surface, i.e. 
in the interval —20 € E < 10 MeV are presented in Fig. 39 for all the Zr-isotopes between A = 108 and 
A — 140. The chemical potential is indicated by a vertical line. For the mass numbers A < 122 the chemical 
potential lies several MeV below the continuum limit (E — 0) and there is only very little occupation in the 
continuum (F > 0). The nucleus ?Zr has a magic neutron number and no pairing. As the neutron number 
goes beyond this closed core, the occupation of the continuum becomes more and more important. Adding up 
the occupation probabilities v? for the levels with E > 0 one finds in the continuum 2 particles for N = 84, 
4 for N = 86,6 N = 88, roughly 3 for N = 90, roughly 4 for N = 92, roughly 3 for N = 94, roughly 4 for 
N = 96, roughly 5 for N = 98, and roughly 6 for N = 100 (where the neutron drip line is reached). 

Considering the very extended neutron rms radii of these systems and estimating the number of valence 
neutrons which would fit into the same volume outside the "??Zr core (shown in Fig. 38 by a dashed dotted 
line) if they would be packed with normal neutron density, one finds for !?^Zr 24 neutrons instead of 12 and for 
1407» 34 neutrons instead of 18. This means that the density outside of the '?Zr core is reduced by roughly a 
factor of 2. This phenomenon is therefore clearly a neutron halo and not a neutron skin. As one sees from Fig. 
3 the tail of the density is proportional to exp(—2vV 4mSr) with S ~ 0.4 MeV. This fact is also emphasized by 
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the extremely low 2n-separation energies (< 0.5 MeV) of these systems. This phenomenon is called a Giant 
Halo because of the large number of particles in the halo region |43].. 


4.4.8 Giant halos in lighter nuclei 


Although the giant halos were predicted in Zr isotopes, it is very difficult for experimentalists to confirm them 
because the exotic Zr isotopes with N > 82 are too heavy to be synthesized by the RIB facilities at present. 
Thus it is valuable to investigate the giant halo phenomena in lighter nuclei which are experimentally accessible 
today. Ground state properties of all the even-even O, Ca, Ni, Zr, Sn, Pb isotopes ranging from the proton 
drip line to the neutron drip line were investigated [44,126]. The two neutron sepation energies are shown in 
Fig. 14 where one finds the Sən values for exotic Ca isotopes near neutron drip line are very close to zero in the 
large mass region, i.e., Son ez 2.06, 1.24, 0.76, 0.43, 0.21, 0.04 MeV for A = 62, 64, 66, 68, 70, 72 isotopes, 
respectively. If one regards 99 Ca as a core, then the several valence neutrons occupying levels above the N = 40 
sub-shell for these A > 60 exotic nuclei are all weakly bound and can be scattered easily to the continuum 
levels due to the pairing interaction, especially for 90772?Ca. This case is very similar to Sən in Zr isotopes 
with N > 82 (see Fig. 14) [43] . Note that for Sn isotopes in the vicinity of the drip line [124], however, Sa, 
decrease rather fast with the mass number and for Ni isotopes [42], Sən are quite large (~ 2MeV) and undergo 
a sudden drop to negative value at the drip line nucleus. Such behavior of S5, in Ca isotopes with N > 40 
indicates that there exist giant halos in these nuclei, just as what happens in Zr isotopes with N > 82 [43]. 
From Fig. 18, one also finds similarity between Ca and Zr isotopes: the neutron radius increases suddenly with 


N neutron number increasing in the Ca chain above N = 40 and in Zr chain above N = 82. Together with the 
— analysis of nucleon density distributions, single particle levels and the occupation probabilities in the canonical 
> ` basis, it is predicted that giant halo may exist in Ca isotopes with A>60 [44, 126]. 


Compared with giant halos in Zr chain, the halos in Ca isotopes will be much easier to access with today's 
radioactive facilities. This means that the giant halos in the Ca isotopic chain are challenging but promising 
. for experimentalists. In order to make the conclusion more solid and exclude the effects from possible error in 


dí Sən (~ 1 MeV), in Fig. 40, the two-neutron separation energies Sən in drip line region for even-neutron Ne, Na, 


Mg, and Al chains (left panel) and for Ar, K, Ca, Sc and Ti chains (right panel) are presented. Open symbols 


c represent the values calculated from the RCHB theory with the NLSH parameter set while corresponding filled 
5. symbols represent the data available. It should be noted that the calculations are performed by keeping the 


æ shape spherical. Nevertheless it can be seen from the figure that all these 55, curves are almost parallel to 


„Z each other and the zero Sən line. Although the calculated Sən value may differ from the data by 1 or 2 MeV, 


(© there must be one or more isotope chains in which many nuclei have the $o„ values very close to zero thus are 


© good candidates for giant halos. 

= Experimentally much efforts have been made to search the drip line for heavier elements. New experiments 
) have shown that both 37Na and *4Ne are bound while ??Ne and ?9Na are missing [190,191], as expected from 

the predictions in Ref. [70]. In Fig. 40, "Na and *4Ne lie in the drip line region and approach the predicted 

giant halo area. Therefore much more experimental effort should be devoted to measuring their masses and 

extending the isotope chains for the possible giant halos. 


4.5 Halos in hypernuclei 


The production mechanisms, spectroscopy, and decay modes of hypernuclear states have been the subject 
of many theoretical studies. The relativistic Hartree Bogoliubov model in coordinate space with finite range 
pairing interaction [71,121] has also been applied to describe A hypernuclei with a large neutron excess [192]. It 
is found that the inclusion of the A hyperon does not produce excessive changes in bulk properties but shifts the 
neutron drip by stabilizing an otherwise unbound core nucleus at drip line. For example, the nucleus *?Ne was 
unbound without the A (see Fig. 41 (a)). The presence of the strange baryon stabilizes the otherwise unbound 
core and ^?*^Ne becomes bound (see Fig. 41 (b)) The microscopic mechanism through which additional 
neutrons are bound to the core originates from the increase in magnitude of the spin-orbit term in presence of 
the A particle [192]. 

Motivated by experimental knowledge of A-N interaction and theoretical understanding on giant halo [43, 
44, 126], it is interesting to observe the possible appearances of halos in hyper exotic nuclei as what have been 
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Figure 40: Left: The two-neutron separation energies Sən for even-neutron Ne, Na, Mg, and Al nuclei in 
the drip line region. Open symbols represent the values calculated with the RCHB theory with the NLSH 
parameter set, while corresponding solid symbols represent the data available. The horizontal line at 2MeV 
denotes the upper limit for the possible occurrence of halos. Right: for even Ar, K, Ca, Sc, and Ti nuclei in 
the drip line region. Taken from Ref. [126]. 
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( J Figure 41: 1f-2p single-particle neutron levels in the canonical basis for the Ne (a), and Ne--A (b) isotopes. 
. The dotted line denotes the Fermi level. Taken from Ref. [192]. 


( ` done in Ref. [47]. In Fig. 42, two neutron separation energies S», for ordinary nuclei, single-A hyper nuclei 
== and double-A hyper nuclei of Ca isotopes labelled by A = 0, A = 1 and A = 2, respectively, from the proton 


£ = to neutron drip line are presented. It is clear that one or two A hyperons lower the neutron Fermi level a little 


' bit but keep the neutron shell structure unchanged. Therefore, the neutron drip line is pushed outside from 
N = 52 in ordinary isotope chain to N = 54 in hyper isotope chain. Meanwhile, giant halo due to pairing 
correlation and the contribution from the continuum still exist in Ca hyper nuclei similar as that in ordinary 
Ca isotopes [44,126]. This is a slight but rewarding step for exploring the limit of drip line nuclei on the basis 
of the giant halo. For more details, see Ref. [47]. 


Apart from neutron halos in hyper nuclei, as A hyperon is less bound than the corresponding nucleon in 
nuclei, it is worth investigating the existence of hyperon halos. In Ref. [194], a lambda halo was tentatively 
suggested in a RMF plus BCS model. However, the normal BCS method will have unphysical solution by 
involving baryon gas [42], such prediction needs further study. Therefore hyper carbon isotopes are studied 
by RCHB theory in Ref. [193]. In Fig. 43, the baryon density distributions in logarithmic scales versus the 
radius r in hyper nuclei SO, zx and DC are presented, where the solid, dot-dashed, long-dashed and dotted 
lines represent density distributions for the lambda, neutron, proton, and corresponding core, respectively. It 
can be seen that up to two A hyperons added to the core !?C, the nucleon density distributions remain the 
same and hyperon density distributions at the tail are comparable with those of the nucleons. An intriguing 
phenomenon appears in AO where the hyperon density distribution has a long tail extended far outside of its 
core 44C (denoted by dotted line), as shown in the right panel of Fig. 43. This is a signature of hyperon halo. 
It is due to the weakly bound state Im in no which has a small hyperon separation energy and a density 


57 


10 20 30 40 50 60 


dripline nuclei 


S,, [MeV] 


S,, [MeV] 


S,, [MeV] 


Figure 42: The two-neutron separation energies Dän in even-N Ca isotopes versus the neutron number N. The 
upper panel is for ordinary nuclei, the middle one for single-A hyper nuclei, and the lower one for double-A 
hyper nuclei. Taken from Ref. [47]. 
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Figure 43: Density distributions for A (solid), neutron (dot-dashed), proton (long-dashed) and corresponding 
core (dotted) in IPC, 54C and 13 C in logarithmic scales. Taken from Ref. [193]. 


distribution with long tail. 


5 Magic numbers for superheavy nuclei 


One of the fundamental and persistent questions in nuclear science is the exploration for the limits of charge and 
mass that a nucleus can attain and the creation of nuclei with masses and charges much larger than those we 
are familiar with, i.e., superheavy nuclei. The properties of superheavy nuclei are important to understand not 
only the nuclear structure, but also the structure of the stars and the evolution of the universe. However, after 
lots of efforts for both higher luminosities and better detection efficiencies as well as the impressive progresses 
on the synthesis of superheavy nuclei, the borders of the upper-right corner in the nuclear chart still remain 
unknown. 

Due to fission resulting from the balance of the attractive nuclear surface tension and the repulsive Coulomb 
force, the superheavy nuclei cannot exist in the classical charged liquid-drop model. The stability of superheavy 
nuclei is mainly determined by shell effects. At the magic proton or neutron numbers 2, 8, 20, 28, 50, and 
82, as well as N=126 for neutrons, nuclei have higher stability and abundance compared with their neighbors. 
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Particularly, the highest stability is observed in the case of the doubly magic nuclei. Amongst other special 
properties, the doubly magic nuclei are spherical. A semblance of the same would work also for superheavy 
nuclei, if there were magic numbers in this region. Consequently, these nuclei will be guarded against a faster 
decay by fission and have more opportunities to be bound. Therefore, it is important to find out the regions 
in the (Z, N) plane where the shell effects are strong and the long lifetimes of the superheavy nuclei can be 
expected. 

The pioneering works on the superheavy elements are a series of macro-microscopic nuclear shell model cal- 
culations performed over more than 30 years to determine nuclear masses and energy surfaces, which suggested 
the existence of a group of relatively stable nuclei separated in neutron and proton numbers from the known 
heavy elements by a higher instability when approaching the closed spherical shells Z=114 and N=184 [195-198]. 
'This group of nuclei become to be known as the island of superheavy elements and the verification of their 
existence is one of the most challenging topics in world-wide heavy ion research facility. 

Experimentally, however, although the superheavy elements up to Z=116 are synthesized or claimed to be 
synthesized [199-212] by different heavy ion reaction types including the cold, warm and hot fusion reactions, 
they are not examples of the originally thought island of superheavy elements yet. 

There is no consensus among different theories with regard to the center of the island of superheavy nuclei. 
Based upon phenomenological models such as finite-range droplet model (FRDM), the shell closures were 
predicted at Z=114 and N=184 [213]. Additionally, the FRDM also predicts larger shell gaps at Z=104, 
106, 108, 110 and at N=162, 164 [213]. In Nilsson-Strutinsky scheme, a similar pattern of deformed nuclei 
have been predicted about Z=108 and N=162 as in FRDM [214,215]. However, the main obstacle is the 


(O question whether the macroscopic approaches which apply to the region of 8-stability line can be extrapolated 
ON to the superheavy nuclei. Recently, microscopic calculations [87, 216-226] are also attempted in describing 
J the superheavy nuclei. In the framework of RHB theory, calculations with a finite-range pairing force of 


€ ^ Gogny interaction D1 and effective interaction NLSH show that Z=114 and N=160, N=166, N=184 exhibit 


JO stability compared to their neighbors and indications for a doubly magic character at Z=106 and N=160 are 


> also observed [217]. The Skyrme Hartree-Fock (SHF) method with interactions SkP and SLy7 predicts magic 


^ numbers at Z=126 and N =184, and also predicts the increased stability due to the deformed shell effects at 


ey N = 162 [218]. Considering non-relativistic SHF effective interactions SkM*, SkP, SLy6, SkI1, SkI3, SkI4 
^ and relativistic mean field (RMF) effective interactions PL-40, NLSH, NL-Z, TM1, the doubly magic spherical 


— nuclei 184114, 172120 and 184126 based on two-nucleon gaps ô2p and ó», are given [219]. The uncertainty on 
“== the magic numbers lies in the uncertain strength of the spin-orbit coupling in non-relativistic approaches or 
^^ the effective interactions in relativistic models. The prediction of magic numbers in superheavy nuclei remains 
; ` a challenge for the nuclear models, which have to be rigorously tested by a wide variety of nuclear properties 


— throughout the periodic table. 

: As the RCHB formalism allows for the proper description of the coupling between the bound states and the 
' continuum by the pairing force and has shown a remarkable success in the description of nuclei with unusual 
N/Z ratio [36, 42-44, 46], it is suitable to apply the RCHB theory for superheavy nuclei. In principle, only a 
calculation in a large multidimensional deformation space can definitively decide the appropriate ground-state 
shape. However, as the traditional superheavy nuclei are expected to be spherical and are located on the 
nuclear chart around a spherical doubly magic nucleus next to DE Pb, the RCHB theory with the assumption 
of spherical shape can be applied for the search. Therefore here we will mainly focus on the magic numbers 
in spherical superheavy nuclei. Some crucial observables in searching shell closures and magic numbers are 
taken into account, including two-nucleon separation energies, two-nucleon gaps, the shell correction energies, 
as well as the pairing energies and the effective pairing gaps. At the end of this section, a brief discussion on 
the stability of the doubly spherical nuclei against deformation and the main progress on the investigation of 
the deformed superheavy nuclei will be given, as well as the a— decay half-lives of superheavy nuclei. 


5.1 Two nucleon separation energies 


Traditionally, the magic nuclei are characterized by higher stability than their neighboring ones and the magic 
numbers can be revealed by the differences of the binding energy, i.e. the one- or two-nucleon separation 
energies [16]. Therefore, the nucleon separation energies are good starting points to investigate the magic 
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numbers in superheavy region. Due to the absence of odd-even effects, the two-nucleon separation energy Sop 
(Sən) is better to quantify shell effects than the single-nucleon separation energy Sp (Sn). For an isotonic 
(isotopic) chain, the $5, (S2n) decrease with the proton (neutron) number and its sudden jump indicates the 
occurrence of a proton (neutron) magic number. 
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Figure 44: The two-proton(left panel) and two-neutron (right panel) separation energies, So, and Son, as a 
function of mass number A obtained by RCHB calculation with PK1 effective interaction. Taken from Ref. [48]. 


d In Ref. [48], the two-proton and two-neutron separation energies have been investigated in the superheavy 
region with proton number Z ranging from 100 to 140 and neutron number N from (Z + 30) to (2Z + 32) 
» by the RCHB theory with NL1, NL3, NLSH, TM1, TW99, DD-ME1, PK1 and PKIR. As an example, the 
results with PK1 are shown as a function of mass number A in Fig. 44. For the proton in left panel, the large 
gaps appear at Z = 106, 120,132,138, while the usual Z=114 does not appear. Furthermore, all the gaps 
are strongly neutron dependent. On the neutron side in right panel, the relatively larger gaps can be seen at 
IN —138, 172, 184, 198, 228, 238, 252, 258, and 274, and also the size and the shape of the gaps differ notably. 
From similar investigations with other effective interactions [48], it is found that the magic proton numbers 
Z=120, 132, and 138 are common while Z=106 is observed only for NL3, NLSH, TW99, PK1, and PKIR. 
* Furthermore the gaps do not exist for all isotopes. On the neutron side, Sən show distinguishable gaps at 
^ N—138, 172, 184, 198, 228, 238, 252, 258, and 274 with all the effective interactions. Apart from these, there 
are also other gaps which appear only for some effective interactions, e.g., N = 164 for NLSH, TW99, DD-ME1, 
PK1 and PKIR, and N = 216 for NLSH, TW99, PK1, and PKIR. For the magic numbers marked either by 
the common gaps or the interaction-dependent gaps, the shell quenching phenomena, i.e., the gaps at N=184, 
198 (NL1), 216 (NLSH, TW99, PK1, and PK1R), 228 (NL1, NL3, NLSH, TM1, PK1, and PK1R), 238, 252, 
and 258 (NLSH, TM1, PK1, and PK1R) appear only for certain Z, are observed. 


5.2 Two-nucleon shell gaps 


'The changes of the two-nucleon separation energies can be quantified by the second difference of the binding 
energies, i.e., the two-nucleon shell gaps: 


óx(N,Z) = So(N,Z) — Sop(N, Z + 2) = 2EB(N, Z) — Ep(N, Z 4-2) — Ep(N, Z — 2) (104) 
5an(N, Z) Son(N, Z) — Son(N + 2, Z) =2Ep(N,Z) — EB(N + 2, Z) — EB(N — 2, Z). 


A pronounced peak of the two-nucleon shell gaps corresponds to a drastic change of the two-nucleon separation 
energies and indicates the shell closure [227]. So far, the two-nucleon shell gaps have been extensively used to 
be an indicator for the magic number [219,221] and to analyze the shell quenching phenomenon [228]. 

The two-proton gaps 02, from the RCHB calculations for even-even nuclei with Z—102 - 138 as a function 
of Z are shown in Fig. 45. The two-proton gaps 62, with the same N are connected as a curve. A peak at 
certain Z in the curve suggests the existence of magic proton number. The sharpness of the peaks represent the 
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Figure 45: The two-proton shell gaps, ó»,, as a function of proton number obtained by RCHB calculation with 
effective interactions NL1, NL3, NLSH, TM1, TW-99, DD-ME1, PK1, and PK1R, respectively. Taken from 
Ref. [48]. 
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Figure 46: The two-neutron shell gaps, den, as a function of neutron number obtained by RCHB calculation 
with effective interactions NL1, NL3, NLSH, TM1, TW-99, DD-ME1, PK1, and PK1R, respectively. Taken 
from Ref. [48]. 
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goodness of the magic numbers while the quenching effects are associated with the bundle of the curves at the 
certain Z. The size of the gaps of two-proton separation energies 55, in Fig. 45 correspond to the magnitude 
of the peaks of two-proton gaps 09, in Fig. 45. It is observed that common magic proton numbers Z=120, 132, 
and 138 exist for all the effective interactions while Z=106 is observed only for NL3, NLSH, TW99, DD-MEI, 
PK1, and PK1R. Furthermore, the peak at Z=114 for NLSH and TW99 and the peak at Z=126 for NL1 are 
also observed, though they are not so obvious as that at Z=120. 

Similar figure for two-neutron gaps óo, are shown in Fig. 46. There are interaction-independent peaks of 
two-neutron gaps do, at N=138, 172, 184, 198, 228, 238, 258, and 274 which are consistent with the conclusion 
from the two-neutron separation energies $5, . Moreover, small peaks are also observed for N=154 (NLSH 
and TW99), N=164 (NL3, NLSH, TW99, DD-MEI, PK1, and PK1R), and N=216 (NLSH, TW99, PK1, and 
PKIR). 


5.3 Shell correction energies 


'The shell correction energy, representing the overall behavior of the single-particle spectra, is another suitable 
observable to identify the shell closure. This quantity, which is derived from the single-particle spectra by 
the Strutinsky procedure [229], is defined as the difference between the total single-particle energy E and the 


smooth single-particle energy E: 
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= d where N(Z) is the particle number, e; is the single-particle energy, A is the smoothed Fermi level determined 
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e E polynomial LA" (z?). The shell correction energy provides an indicator about the deviation in the level structure 
= of nuclei away from uniformly distributed ones. A large negative shell correction energy corresponds to shell 


>= closure at certain nucleon number. 

i With the proton single-particle spectra in canonical basis from the RCHB calculations, the shell correction 
energies EP, ,1 for even-even nuclei with Z ranging from 100 to 140 are shown as a function of Z in Fig. 47. The 
nuclei with the same N are connected as a curve. A deep valley at certain Z on a curve hints the shell closure. 
The valleys at Z=120, 132, and 138 are conspicuous for all the effective interactions. Possible shell closure 
at Z=106 for NL1, NL3, NLSH, TW99, DD-MEI, PK1, and PKIR, at Z=114 for NLSH, TW99, PK1, and 
PKIR, and at Z=126 for NL1 are also observed. The shell closure from the shell correction energies for proton 
B o are consistent with the peaks of the two-proton gaps 49,, though the magnitude of magicity slightly 
differs. Similar to the two-proton gaps 0», the spread of the valleys of the shell correction energies for proton 
FP, ey indicates the quenching of the magic number. In view of the shell correction energies for proton Ee, ell 
the quenching phenomena of magic proton numbers in RCHB calculation are universal. These facts suggest 
that the magic proton numbers depend also on neutron number N. 

The corresponding shell correction energies for neutron E% .,, for even-even nuclei with N=130 - 312 are 
demonstrated as a function of N in Fig. 48. As the neutron numbers extend from N=130 to N=290, the 
valleys for shell correction energies for neutron EN Am Fig. 48 are not so obvious as those for shell correction 
energies for proton EP, in Fig. 47. However, similar conclusions as the two-neutron separation energies Son 
and the two-neutron gaps da, for magic neutron numbers can also be drawn. Similar calculation has also been 
done in Ref. [230] with SHF and RMF. 
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Figure 47: The shell correction energies for proton, EP, a, as a function of proton number obtained by RCHB 
calculation with effective interactions NL1, NL3, NLSH, TM1, TW-99, DD-ME1, PK1 and PKIR, respectively. 
Taken from Ref. [48]. 
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Figure 48: The shell correction energies for neutron, E% eg as a function of neutron number obtained by RCHB 


calculation with effective interactions NL1, NL3, NLSH, TM1, TW-99, DD-ME1, PK1 and PK1R, respectively. 
Taken from Ref. [48]. 
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The first valley locates in N=184 for all the effective interactions and involves N=172 (NL3, NLSH, TM1, 
TW99, DD-ME1, PK1, and PK1R), even N=164 (NLSH and TW99). And the other valley is the mixture 
of the valleys at N=228 (TW99 and DD-ME1), N=252 (NLSH, TM1, PK1, and PK1R) and N=258. These 
fine structures suggest the smearing of magic neutron numbers. Based on the shell correction energies, the 
shell closures are smeared compared with those from the two-nucleon separation energies 55, and 55, and the 
two-nucleon gaps 9o, and ó»,. However, it is found in Ref. [48] that the magic numbers at Z=120, 132, and 
138 are common but strongly quenched, the shell closures near N=172 and N=184, and N=228, N=252, and 
IN —258 are blurred. The other magic numbers from the shell correction energies, such as Z=106, 114, and 126, 
appear only for some effective interactions. 

In the usual Hartree or HF mean field level, as the nonphysical particle continuum are involved, much 
efforts have been made to avoid the divergence of single particle energy density around the threshold [231-234]. 
Based on the Green's-function approach to the level density, a method of calculating shell corrections has 
been adopted within the SHF and RMF theories to a large-scale survey of spherical shell energies throughout 
the whole landscape of conceivable superheavy nuclei [230]. From SHF with SkP, SLy6, SkI3 and SkI4 and 
RMF with NL3 and NL-Z2, the total shell correction energies Eshe for spherical superheavy nuclei around 
the expected island of stability around Z = 120 and N = 180 are given in Fig. 49. Instead of narrow stripes 
of large Espey localized around magic numbers for normal nuclei, a wide area of shell stabilization is found 
to spreads over all shell closures predicted. In other words, the predicted magic shells for superheavy nuclei 


= are much mellowed, which are consistent with the results from RCHB theory. The disappearance of a familiar 


> pattern of magic numbers and the appearance of broad valleys of shell stability are due to the rather large 
. single-particle level density and the appearance of many low 7 shells around the Fermi level [230]. 
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Figure 49: Total shell correction energies, calculated from SHF with SkP, SLy6, SkI3 and SEL and RMF with 
NL3 and NL-Z2, for spherical even-even superheavy nuclei around the expected island of stability around Z — 
120 and N = 180. The thick solid lines denote two-particle drip lines. Black squares mark nuclei calculated 
to be stable with respect to 8 decay. White color indicates nuclei with positive shell corrections, black color 
denotes nuclei with Eshe beyond -12 MeV. Taken from Ref. [235]. 


5.4 Magic numbers and deformations 
5.4.1 Magic numbers for spherical superheavy nuclei 


After the detailed analysis of the two-nucleon separation energies, 55, and San, the two-nucleon shell gaps, 
Õ2p and äu, the shell correction energies, E? ey and Eu, as well as the pairing effects, the magic num- 
bers for superheavy nuclei in the RCHB theory are found to be [48]: Z = 120,132,138 for proton and 
N = 172,184, 198, 228, 238,258 for neutron. As many theoretical methods have been adopted to probe the 
superheavy magic numbers and there is no consensus among them with regard to the center of the island of 
superheavy nuclei. It is interesting to briefly review the up-to-date studies of the superheavy magic numbers 
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including the present one: 


e Z=106: As we have seen, this shell closure is strongly supported by the effective forces NL3, NLSH, 
TW-99, PK1, and PK1R, partially by NL1 and DD-ME1, but not at all by TM1. We also noted that the 
shell closure at Z=106 depends sensitively on the neutron number and is visible only around N ~ 154 
and ~ 218, while the later can only be clearly seen in the results calculated with TW-99 and NLSH. The 
RHB calculations with the effective force NLSH and the Gogny interaction D1 for the pairing correlation 
predicts a doubly magic character at Z=106 and N=160 [217]. In addition, the shell closures at Z=104, 
106, 108, and 110 and N=162, and 164 which correspond to prolate deformations are predicted by the 
FRDM [213] and the deformed shell closure occurs at Z—108 in Nilsson-Strutinsky scheme [214, 215]. 


e Z—114: The shell closures are predicted at Z=114 and N=184 in the phenomenological shell models [197, 
198,213,214]. However, as we have noted, in the RCHB theory the magicity of Z=114 is not at all obvious, 
in agreement with other RMF investigations [22,217,219,236,237]. By adopting a new effective interaction 
NL-RA1, the nucleus 784114 is predicted to be the next spherical doubly magic nucleus in a RMF plus 
BCS calculation [238]. Within the SHF approach, only SkI4 prefers the doubly magic nucleus at Z—114 
and N=184, which is not supported by other Skryme forces such as SkM*, SkP, SkI1, SkI3, SLy6, and 
SLy7 [219, 221]. 


e Z=120: The shell closure at Z—120 is strongly supported in the RCHB calculations. The doubly magic 
characteristic of 172120 has been addressed within the RMF theory [219, 221, 234, 238] as well as the 
effective field theory [237]. Most Skyrme HF investigations also suggest that this nucleus is doubly magic 
except for SkP [219, 221]. Due to the strong shell quenching effect in the Z—120 isotopic chain, other 
doubly magic candidates with N=184, 198, 228, 238, 252, 258, and 274 are interaction-dependent and 
need to be investigated in more detail. 


e Z=126: Several Skyrme HF investigations with the effective interactions such as SkP, SkM*, SLy6 and 
SLy7 predict doubly magic system at Z=126 and N=184 [218,219,221]. However, this conclusion is not 
supported by the RMF calculations although N=184 might be a neutron magic number. 


e Z=132 and 138: Little attention has been paid to the possible proton shell closures at Z > 130 due 
to their large Coulomb repulsive potentials. However, the shell closures at Z—132 with N ~ 210-240, 
and at Z=138 with N ~ 180-200 and N ~ 258 appear in most relativistic and non relativistic mean field 
models [48, 219], . 


l Normally the magic numbers correspond to the large energy gaps between the single particle levels and the 
‘© doubly magic nuclei are those with protons and neutrons filled up to the gap. It is therefore interesting to 
) study the single-particle level structure in these doubly magic nuclei, or vice versa. 
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Figure 50: Spherical single particle spectra in 292120 for protons (left) and neutrons in various mean-field 


calculations as indicated. Taken from Ref. [221]. 


As examples, the spherical single particle spectra in ?9?120 for protons (left) and neutrons in various mean- 
field calculations are given in Fig. 50 [221]. The occurrence of the shell closure at Z — 120 depends on the 
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amplitude of the spin-orbit splitting of the 3p states above the Fermi level and the 2f levels below the Fermi 
energy. It has also been pointed out that the shell closures may be coupled to the density profiles (and vice 
versa) based on the systematic analysis of density distribution of superheavy nuclei. For 797120, a significant 
central depression or central semibubble is obtained in the spherical case [221,239]. Recently, the influence of 
the central depression in the density distribution of spherical superheavy nuclei on the shell structure is studied 
in some detail within the RMF theory [240]. It is found that a large depression leads to the shell gaps at the 
proton Z = 120 and neutron N = 172 numbers, whereas a flatter density distribution favors N = 184 and leads 
to the appearance of a Z = 126 shell gap and to the decrease of the size of the Z = 120 shell gap. On the other 
side, the deformation effect has been found to be significant in the calculation of the density distribution in 
this nucleus, namely only a weak central depression in the deformed case for 797120 compared to the predicted 
semibubble at the spherical shape [240, 241]. 


5.4.2 Stability of the doubly magic superheavy nuclei against deformation 


The magic numbers discussed above are mainly based on the assumption of spherical geometrical configuration. 
— However, such assumption is not always true for superheavy nuclei [213-215,242]. It is worthwhile to investigate 
= the potential energy surfaces to see the validity of spherical configuration. 


In Fig. 51, the potential energy surfaces for 797120, 304120, 315150, 345150, 355120. and 375120 (solid lines) 


( l obtained from constrained RMF calculations with the pairing correlations treated by BCS approximations 


- are displayed. Although these nuclei have the same proton number, their potential energy surfaces are quite 
different from each other as shown in Fig. 51. For nuclei 797120, 304120 and 378120, there is an obvious local 
minimum with the spherical configuration £5 ~ 0, while another local minimum with large deformation £5 ~ 0.6 


N - can be also clearly seen. For ?/5120, the spherical local minimum is very shallow, and for 348120 and 358120, 


the spherical minimum is hardly seen and a local minimum with f$ ~ 0.25 appears instead. There are also 


© local minimum around 62 ~ —0.4 with oblate deformation. In addition to the local minima discussed above, 
““ we would like to make a few remarks on the absolute minima for these nuclei. For 292120 and 378120, the two 
“> minima with different deformations have almost the same energies, i.e. the so-called “shape coexistence” may 
SS exist. In particular, the spherical minimum is indeed the absolute minimum for 797120 (Eg = -2064.3 MeV 


e with 82 ~ 0 vs. Eg = -2063.6 MeV with 82 ~ 0.6), while the ground state of 378120 is quite delicate with Eg = 
( e -2398.7 MeV for bz ~ 0 and Eg = -2398.7 MeV for Bo ~ 0.6. For 24120, the absolute minimum at Ba ~ 0.6 
= is much deeper (~ 5.5 MeV) than the spherical configuration. In addition for 348120, the absolute minimum 


F ; lies at 62 ~ 0.2, and it is a well-deformed nucleus according to the present calculation. From the RCHB and 


the constrained RMF+BCS calculations, it is shown that spherical doubly magic nuclei may exist, e.g., 797120 
and 378120 [48]. 


The role of shell effects in superheavy nuclei can be also seen in Fig. 51, in which the corresponding 
macroscopic energy, defined as the difference between the binding energy and the total shell correction energy, 
Ep — Eshei, is plotted with dashed lines. It can be seen that without the shell effects, the superheavy nuclei 
hardly exist. As noted long time ago, the shell effects play an essential role to stabilize the superheavy nuclei 
against the fission. 


Apart from axial symmetric case, the influence of other deformation such as triaxiality and reflection 
asymmetry deformations has been investigated as well, e.g., see Ref. [23,226]. Figure 52 presents an example 
of potential energy surface of the typical actinide nucleus ?!?Pu, which has often served as a benchmark for 
mean-field models, obtained with the Skyrme force SLy6. A typical doubled-humped fission barrier is well 
shown. The inner barrier explores triaxial degrees of freedom, which reduce the axial barrier by about 3 MeV, 
while the outer barrier explores reflection asymmetric shape [226]. The similar behavior is held for heavier 
nuclear system [220]. It is also noted that though the further deformation degrees of freedom reduce the fission 
barriers obviously, they hardly change the energies of mimina. Therefore, the constraint RMF model with axial 
symmetry remains valid to confirm the spherical configuration of a superheavy nucleus. 
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Figure 52: Double-humped fission barrier of the typical actinide nucleus Pu. The dotted line denotes 
an axial and reflection-symmetric calculation, the full line denotes a triaxial (inner barrier) and axial and 
reflection-asymmetric calculation (outer barrier). The various shapes along the axial paths are indicated by 
the contours of the total density at ro = 0.07 fm. Taken from Ref. [226] 
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5.5 Alpha decay and the half-lives 


The preferred decay mode of shell stabilized superheavy nuclei is a decay. A key quantity for o decay, the Qa 
value is defined as 


Qal(N, Z) = HUN, Z) - E(N — 2, Z — 2) — E(2,2), (107) 


which is the energy of the a-particles emitted by radioactive heavy and superheavy nuclei. It is interesting to 
compare the Qa data of the discovered superheavy nuclei with predictions from mean-field models. Such works 
have been done in various of mean field framework, e.g. Refs. [23, 222-225, 243]. In the recent experiments 
designed to synthesize the element 115 in the 747Am+“°Ca reaction at Dubna in Russia, three similar decay 
chains consisting of five consecutive a-decays, and another different decay chain of four consecutive a-decays 
are observed, and the decay properties of these synthesized nuclei are claimed to be consistent with consecutive 
a-decays originating from the parent isotopes of the new element 115, ?55115 and ?57115, respectively [212]. The 
deformed RMF+BCS calculation with a density-independent delta-function interaction in the pairing channel 
and the effective interaction TMA is made to study these newly synthesized superheavy nuclei 785115, ?57115, 
and their a-decay daughter nuclei [225]. The calculated a-decay energies and half-lives agree well with the 
experimental values and with those of the macroscopic-microscopic FRDM+FY and YPE+WS models. 

For an area of enhanced stability, the a-decay half-lives are expected to be longer than its neighbors. The 
half-lives of a-decay can be obtained with Qa value by phenomenological Viola and Seaborg systematics [244]: 


logioTy = (1.66175Z — 8.5166)Q2'/? — (0.20228Z + 33.9069) (108) 


) where Z is the proton number of the parent nucleus, T4 in second and the parameters taken from Ref. [245]. 
\ It should be noted that Eq. (108) is based on the WKB approximation, and provides only a rather crude 
` estimation of Ty since it disregards many structure effects such as deformation, and configuration changes, etc. 
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Figure 53: The a-decay half-lives Tą as a function of neutron number obtained by RCHB calculation with PK1 
effective interaction. Taken from Ref. [48]. 


In Fig. 53, the half-lives 74 from the RCHB calculation in logarithm scale are plotted as a function of 
neutron number N with the effective interaction PK1. The half-lives corresponding to 1 ns, 1 us, 1 ms, 1 s, 
1 h, 1 y, and 1 ky are marked by the dashed lines. Each curve in this figure corresponds to an isotopic chain 
in the region with Z ranging from 102 to 140. In Fig. 53, the half-lives logıo7„ increase with N. The jumps 
of the curves correspond to magic proton numbers Z=106, 120, 132, and 138, which depend on the effective 
interactions. The peaks of each curve correspond to magic neutron numbers, i.e., N —172, 184, 228, 238 and 
258. It can be seen that the magic numbers suggested by two-neutron separation energies S5, (two-neutron gaps 
dyn) or two-proton separation energies $5, (two-proton gaps ô2p) can also be found here. It should be noted 
that the half-lives Tẹ shown here are obtained from the spherical RCHB calcualtions. In fact, both the pairing 
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correlation and deformation are essential to reproduce the experiment Qa values [222,224,246]. Calculations 
along this line are necessary. 


6 Summary and perspectives 


This article has been concerned with the relativistic description of the exotic nuclei particularly with the 
Relativistic Continuum Hartree-Bogoliubov (RCHB) theory and its applications. The last two decades have 
witnessed the development of the relativistic mean field (RMF) theory and its successful application in nuclear 
physics. At the very beginning, the question on the clear signal of relativistic effects in atomic nuclei has puzzled 
us for quite a long time. Now after lots of efforts, it become evident that there are at least the following points 
which do favor the relativistic description of nuclei: 


e One may say that the kinetic energy of a nucleon moving in an nucleus is much smaller than its mass. 
However, as far as the mass is concerned, we have to remember that the mass of a nucleus is smaller 
than the sum of the mass of its consistent nucleons, which is naturally obtained in a relativistic model. 
Furthermore there are always a large scalar and a large vector potentials which are both comparable to 
the mass of the nucleons in nucleus. 


e Relativistically the spin-orbit potential comes out naturally. Furthermore the long existing problems on 
the origin of pseudo-spin symmetry were solved. A very well developed spin symmetry in single anti- 
neutron and single anti-proton spectra is found. This spin symmetry in anti-particle spectra and the 
pseudo-spin symmetry in particle spectra have the same origin. 


e In the framework of the relativistic approach, the nucleons interact via the exchanges of mesons and 
photons. From this point of view, the relativistic approach is more microscopic and fundamental in the 
sense of describing the nuclear system at the meson level. 


e In practice, it turns out that the RMF theory can reproduce better the nuclear saturation properties (the 
Coester line) in nuclear matter and the isotopic shifts in the Pb-region, and is more reliable for nuclei far 
away from the line of 6-stability. 


== Altogether the relativistic mean field theory is a reliable and effective model for nuclear structure. Therefore, 


‘= with the fast progress in producing the exotic nuclei, it is natural to describe the exotic nuclei based on the 

= RMF theory. 

es The formalism, the numerical solutions and the effective interactions for the RMF theory as well as its 

` > application for nuclear matter and neutron stars are presented in Section 2. The numerical solutions of the 
- Dirac equation for spherical nuclei, particularly on the Woods-Saxon basis, are discussed in detail. The new 

effective interactions with nonlinear self-coupling meson fields and density dependent meson-nucleon couplings, 

and their influences on properties of nuclear matter and neutron stars are also discussed. 

Pairing correlation plays important roles in open shell nuclei. The conventional BCS method can be 
combined easily with the RMF theory, but it is not justified for exotic nuclei because it could not include 
properly the contribution of continuum states. The formalism and related discussions of the RCHB theory, in 
which the coupling between the bound states and continuum could be included self-consistently, are given in 
Section 3. Discussion on the continuum states, the resonant BCS method and several methods for obtaining 
single particle resonant states are briefly reviewed. 

The study of exotic nuclei has attracted world wide attention due to their large N/Z ratios (isospin) and 
interesting properties such as halo and skin. In Section 4 we reviewed the application of the RCHB theory to 
exotic nuclei, the binding energies, particle separation energies, the radii and cross sections, the single particle 
levels, shell structure, the restoration of the pseudo-spin symmetry, the halo and giant halo, and halos in hyper 
nuclei, etc. For examples, (1) the Li isotopes were investigated and the halo phenomenon in !!Li was found 
to be the result of the scattering of particle pairs into continua; (2) the giant halos were predicted in Zr, Ca 
and even lighter nuclei when neutron drip line is approached based on the analysis of two neutron separation 
energies, single particle energy levels, the orbital occupation, the contribution of continuum and nucleon density 
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distribution. These predictions are supported to some extent by the recent experimental information; (3) exotic 
phenomena in hyper-nuclei were also studied; etc. 

The RCHB theory has been applied to study superheavy nuclei which is reviewed in Section 5. The possible 
doubly magic superheavy nuclei have been searched within the RCHB theory using a variant of the effective 
interactions. Base on the detailed analysis on the two-nucleon separation energies 55, and 5»,, two-nucleon 
shell gaps do, and ôn, the shell correction energies E^, ;, and Eu: as well as the pairing energies F air and 
Ei, the effective pairing gaps A, and An, the proton and neutron shell closures have been predicted by using 
the RCHB theory with the effective interactions NL1, NL3, NLSH, TM1, TW-99, DD-ME1, PK1 and PKIR. 
Proton numbers Z=120, 132, and 138 and neutron numbers N=172, 184, 198, 228, 238, and 258 are supported 
by all effective interactions to be magic. The deformation-constrained RMF calculations have been done for 
these doubly magic nuclei candidates and discussions on the deformation of these nuclei are also presented. 

Among many perspectives of the application of the RMF and the RCHB theories, we mention briefly the 
following: 


1. Attempts have been made by extending the present RCHB theory for the axially deformed nuclei. Al- 
though the coupled Dirac equations and Hartree-Bogoliubov equations with fixed pairing potentials have 
been solved in coordinate space, the simultaneous self-consistent solution of axially deformed RCHB equa- 
tions is proved to be very difficult [247]. Therefore the Woods-Saxon basis has been suggested to solving 
the RMF theory in order to generalize it to study exotic nuclei by solving either the Schródinger equation 
or the Dirac equation [73]. It combines the advantages of the efficiencies in harmonic oscillator basis and 
the appropriateness for exotic nuclei in coordinate space. For spherical case, the Woods-Saxon basis has 
been proved to be very successful [73]. The development of deformed RCHB in the Woods-Saxon basis 
is in progress. 


2. The Shell-model-Like APproach (SLAP) method has been suggested to treat pairing correlation in RMF 
theory [248]. With the single particle bound states from the RMF calculation and resonant states obtained 
from RMF theory combined with the Analytic Continuation of Coupling Constant (ACCC) method [110, 
111] or the phase shift scattering method [41,117], the SLAP for the pairing can also be used to describe 
exotic nuclei as well. 


3. Since the pion does not contribute in the Hartree level, we did not mention its contribution throughout 
the paper. However, if the Fock term is considered, one must include pion in the Lagrangian. This would 
bring out some complication in solving the equations as well as new and interesting physics. 


4. Many symmetries have been broken in the mean field theory. Restoration of some symmetries could 
be achieved by making projections or corrections. The microscopic center of mass correction for the 
broken translation invariance in the RMF model has been discussed [51]. The single particle states of a 
deformed mean field do not have good total angular momentum, so is the total intrinsic wave function 
of the deformed nucleus. Thus one must project out a state with good total angular momentum from 
the intrinsic wave function. This can be done by the full angular momentum projection method [35] or 
approximately in a simple way [249]. Recent works on the restoration of rotation symmetry |250, 251] 
are very interesting. Possible extensions and applications for chiral doublet bands and nuclear isomeric 
states are expected [252-258] . 


5. In most of the RMF calculations, the no-sea approximation is applied. Therefore the anti-nucleon degrees 
of freedom are not considered. It would be quite interesting to investigate these degrees of freedom and 
the polarization of the vacuum as well as abnormal nuclei containing anti-baryons [33, 47, 193, 259, 260]. 
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